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We present an easy-to-use analytic toolbox for the calculation of short-time transport properties of 
concentrated suspensions of spherical colloidal particles with internal hydrodynamic structure, and 
direct interactions described by a hard-core or soft Hertz pair potential. The considered dynamic 
properties include self-diffusion and sedimentation coefficients, the wavenumber-dependent diffusion 
function determined in dynamic scattering experiments, and the high-frequency shear viscosity. The 
toolbox is based on the hydrodynamic radius model (HRM) wherein the internal particle structure is 
mapped on a hydrodynamic radius parameter for unchanged direct interactions, and on an existing 
simulation data base for solvent-permeable and spherical annulus particles. Useful scaling relations 
for the diffusion function and self-diffusion coefficient, known to be valid for hard-core interaction, 
are shown to apply also for soft pair potentials. We further discuss extensions of the toolbox to 
long-time transport properties including the low-shear zero-frequency viscosity and the long-time 
self-diffusion coefficient. The versatility of the toolbox is demonstrated by the analysis of a previous 
light scattering study of suspensions of non-ionic PNiPAM microgels [Eckert et al, J. Chem. Phys., 
2008, 129 , 124902] in which a detailed theoretical analysis of the dynamic data was left as an open 
task. By the comparison with Hertz potential based calculations, we show that the experimental data 
are consistently and accurately described using the Verlet-Weis corrected Percus-Yevick structure 
factor as input, and for a solvent penetration length equal to three percent of the excluded volume 
radius. This small solvent permeability of the microgel particles has a significant dynamic effect at 
larger concentrations. 


I. INTRODUCTION 

Suspensions of globular colloidal particles with internal 
hydrodynamic structure, and different surface boundary 
conditions, are abundant in soft matter science. Exam¬ 
ples of technological and biomedical relevance are non¬ 
ionic and ionic microgel particles, and core-shell particles 
consisting of a dry spherical core and a shell of some soft 
material such as a polymer brush. These particulate sys¬ 
tems are to a certain degree permeable to the solvent. Mi¬ 
crogel particles in particular consist of a network formed 
by cross-linked polymer chains. They have useful features 
such as temperature-, pH-, salinity- and concentration- 
dependent [lj swelling behavior as well as elasticity and 
flexibility. This renders them as good candidates for var¬ 
ious applications such as drug delivery agents M , the 
engineering of tissues @-@|, and the modification of rhe¬ 
ological properties [9:]. The elasticity of microgels can 
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be controlled, e.g., by the amount of crosslinker, and the 
length of polymer chains used in the synthetization pro¬ 
cess. Microgels can be therefore considered as bridging 
the gap between genuine hard spheres and ultra-soft col¬ 
loids flpj ]. 

Although microgel and core-shell particle systems have 
been intensely studied experimentally over the past years, 
a quantitative theoretical description of the diffusion 
and rheological properties of concentrated suspensions 
is still on demand. This is owed to the complicated 
many-particle hydrodynamic interactions (His) which 
are significantly influenced by the hydrodynamic struc¬ 
ture of the particles. A theoretical understanding of 
the influence of His on colloidal transport properties 
such as translational and rotational diffusion coefficients, 
the generalized sedimentation coefficient (hydrodynamic 
function), and high-frequency and zero-frequency viscosi¬ 
ties is of key importance also in process engineering, e.g. 
in filtration and fractionation processes [ll|, [l2| , and for 
the energy cost reduction in the transportation of col¬ 
loidal suspensions through viscosity minimization. 

On a coarse-grained level, the porosity-averaged fluid 
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flow inside a solvent-permeable particle is commonly de¬ 
scribed by the Brinkman-Debye-Bueche (BDB) equation 
invoking the Darcy permeability, re 2 , where 1/re is the 
hydrodynamic penetration length BED. Globular par¬ 
ticles with an on average spherically symmetric hydro- 
dynamic structure can be characterized by a permeabil¬ 
ity coefficient, re(d), depending on the radial distance, 
d, from the particle center ITbl . Versatile hydrodynamic 
simulation tools such as the HYDROMULTIPOLE hy¬ 
drodynamic force multipole method m have been de¬ 
veloped which allow for calculating transport properties 
of concentrated dispersions of hydrodynamically struc¬ 
tured particles with the full inclusion of His. However, 
these simulations are numerically expensive, and in prin¬ 
ciple they must be performed separately for each particle 
model. In a recent series of papers, various short-time 
dynami c prope rties of dispersions of uniformly permeable 
spheres [17H20| , and of core-shell particles with uniformly 
permeable shell [2l1423i | . have been calculated using the 
HYDROMULTIPOLE method, as functions of particle 
concentration, reduced Darcy permeability, and shell- 
thickness to particle size ratio. The non-hydrodynamic 
direct particle interactions in these simulations have been 
taken for simplicity as pure hard-core interactions charac¬ 
terized by the excluded volume particle radius a = <r/2. 
Any softness in the effective pair potential, V(r), be¬ 
tween two globular particles at center-to-center distance 
r is hereby disregarded. 

As discussed in mm, a simplifying concept allowing 
for abstracting from specific intraparticle structures is 
the so-called hydrodynamic radius model (HRM) which 
invokes the notion of an apparent no-slip hydrodynamic 
particle radius ah (see also [23, B ). The HRM amounts 
to approximating a globular particle of spherically sym¬ 
metric hydrodynamic structure by a no-slip sphere of 
hydrodynamic radius ah while leaving the effective pair 
potential unchanged. Under from an experimental view¬ 
point surprisingly general conditions, ah is unequivocally 
determined from the measurement of a single-particle 
transport property such as the translational diffusion co¬ 
efficient Dq or the intrinsic viscosity [77] . The definition of 
the HRM includes also spherical particles with fuzzy hy¬ 
drodynamic structure and no sharp outer boundary, and 
with a soft pair potential such as for weakly cross-linked 
ionic microgels jl, B, [H)- For spherical particles hav¬ 
ing excluded volume interactions only with ah < a , the 
HRM reduces to the so-called spherical annulus model. 
For the annulus model, numerically precise simulation 
results for various short-time dynamic properties have 
been given in [2lj . The good accuracy of the simplifying 
HRM was demonstrated in JjLEj, [2lJ, [29( for uniformly per¬ 
meable and core-shell spheres with pure excluded volume 
interactions, by a thorough comparison with simulation 
results. While a single hydrodynamic radius suffices to 
characterize the hydrodynamic intraparticle structure of 
many experimentally realized suspensions, regarding its 
influence on configuration-averaged transport properties, 
the replacement of the soft pair potential by an effective 


hard-core potential is in general a less successful strat¬ 
egy. Methods of calculating static suspension properties 
based on an effective hard-sphere potential such as the 
Barker-Henderson perturbation scheme, a second virial 
coefficient mapping, and additional variational methods 
commonly fail if the longer-ranged, soft part of the pair 
potential stretches out significantly beyond the physical 
excluded volume radius. For c harg e-stabilized colloids, 
e.g., this has been shown in [301 [3l[. 

In this article, we present an easy-to-apply set of ana¬ 
lytic methods, referred to for short as a toolbox, to calcu¬ 
late short-time transport properties of concentrated dis¬ 
persions consisting of spherical colloidal particles with 
internal hydrodynamic structure, and with direct inter¬ 
actions given by the hard-core and two-parameter soft 
Hertz potentials. The latter potential is continuous and 
bounded, and it describes the energy penalty cause d by 
the elastic deformation of two colliding spheres [27], [32] ■ 
The Hertz potential has been shown to be a useful de¬ 
scription of the direct interactions between non-ionic, soft 
microgels of low cross-link density, and this even though 
the particles are of a distinctly inhomogeneous structure 
[23, (H). Our toolbox is based on the HRM, and it takes 
advantage of the tabulated simulation data for spherical 
annulus particles listed in |Tlj]. The toolbox incorporates 
in particular useful approximate scaling relations for the 
wavenumber-dependent sedimentation coefficient, H(q ), 
the short-time self-diffusion coefficient Ds, and the high- 
frequency viscosity These quantities are routinely 
determined in dynamic scattering experiments. The scal¬ 
ing relation expressions are known for permeable parti¬ 
cles with hard-core interactions to be in remarkably good 
agreement with simulation data B> EM H3- We show 
that they apply likewise to particles with a soft pair po¬ 
tential, and we augment them by scaling expressions for 
the collective diffusion coefficient, Dc, and the associated 
sedimentation coefficient K. 

Moreover, we extend the toolbox to long-time dy¬ 
namic properties of concentrated systems of hydrody¬ 
namically structured particles, including the low-shear 
zero-frequency suspension viscosity 77 , and the long-time 
translational self-diffusion coefficient Dl. Different from 
their short-time siblings, long-time transport proper¬ 
ties are affected additionally by the non-instantaneous 
microstructural relaxation of the cloud of neighboring 
Brownian particles. This relaxation is controlled both 
by direct and hydrodynamic interactions. The toolbox 
extension to long-time properties combines the HRM 
with a factorization approximation method introduced 
originally by Medina-Novola [35), and elaborated subse¬ 
quently by Brady [36j, [37] and Banchio et al. 38) . The 
HRM is useful also for long-time properties, since the hy¬ 
drodynamic mobilities in the generalized many-particle 
Smoluchowski diffusion equation describing the configu¬ 
rational distribution function are time-independent [311 ] . 

We demonstrate the accuracy of our user-friendly tool¬ 
box through the analysis of a light scattering study on 
a concentration series of non-ionic, submicron-sized PNi- 
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PAM microgel particles dispersed in dimethylformamide 
(DMF). We show that the static and dynamic scatter¬ 
ing data for the static structure factor S(q) and hydro- 
dynamic function H(q), respectively, can be quantita¬ 
tively described in the complete experimental wavenum¬ 
ber range, on using a deduced solvent penetration length 
equal to three percent of the particle diameter. Calcu¬ 
lations based on the Hertz potential show that the PNi- 
PAM microgels behave statically as effective hard spheres 
in the whole fluid-phase concentration regime. 

The paper is organized as follows: The structure fac¬ 
tor and the associated radial distributon function (RDF), 
g{r), of the Hertz potential model are discussed in Sub¬ 
sec. Ill Al and compared with static light scattering data 
of non-ionic PNiPAM microgels. Subsec. Ill Bl describes 
how hydrodynamic particle structures are related to the 
generic HRM, exemplified for the two important exam¬ 
ples of uniformly permeable spheres and non-permeablc 
rigid spheres with partial hydrodynamic surface slip. 
Sec. EH includes the discussion of various short-time dy¬ 
namic properties, including H(q ), Dg, K, and the high- 
frequency viscosity r] qq. Approximate analytic expres¬ 
sions for these quantities are presented, and the accuracy 
of these expressions is scrutinized against simulation re¬ 
sults for the spherical annulus model. In Sec. Ml the 
toolbox results are compared with static and dynamic 
light scattering (SLS and DLS) data by Eckert et al. on 
PNiPAM microgel suspensions. Sec. [V] describes the ex¬ 
tension of the toolbox to long-time transport properties 
including the long-time translational self-diffusion coeffi¬ 
cient Dl and the zero-frequency viscosity q. Our conclu¬ 
sions are contained in Sec. eh The pairwise additivity 
(PA) and the self-part corrected Beenakker Mazur meth¬ 
ods of calculating short-time transport properties are ex¬ 
plained in the Appendices A and B, respectively, in the 
context of the HRM. 


II. EQUILIBRIUM MICROSTRUCTURE AND 
HYDRODYNAMIC RADIUS MODEL 


A. Static pair correlations 


Owing to the many possible applications, various the¬ 
oretical schemes have been developed for the analytic 
calculation of static properties of microgel suspensions. 
Progress in this direction was made, in particular for 
ionic microgels, through the development of effective pair 
potentials characterized by the suspension temperature 
and salinity, and the bare charge of the microgel particles 
EoHuT ] . The validity of these effective pair potentials has 
been scrutinized in various joint theoretical-experimental 
studies (see e.g. [53, E5J). In contrast to ionic micro- 
gels, the effective pair potentials used for non-ionic soft 
microgels where short-range interactions are not masked 
by the longer-ranged electrostatic repulsion are to date 
still on a more heuristic level. The pronounced depen¬ 
dence of the pair potential in non-ionic microgel systems 


on ambient and intra-particle conditions such as the sol¬ 
vent quality and temperature, number and distribution 
of crosslinker, and length and functionality of polymer 
chains requires a larger number of parameters character¬ 
izing the effective pair potential on a microscopic level. 
On a more coarse-grained level, simplifying pair poten¬ 
tials have been used such as the hard-sphere [43| and 
elastic Hertz HzL EH potentials, and certain ultra-soft 


pair potentials [ 101 ]. The intricate dependence on envi¬ 
ronmental parameters is hidden in these coarse-grained 
potentials in a reduced number of interaction parameters 
such as the effective interaction strength and the effective 
particle radius. 

A useful coarse-grained effective pair potential for non¬ 
ionic globular microgel particles is the Hertz potential, 


PV(r) = 




( 1 ) 


r > a s , 


which describes the elastic deformation of two spheres in 
contact. Here, /3 = l/(fcsT) is the reduced inverse tem¬ 
perature with Boltzmann constant k B and absolute tem¬ 
perature T, and er s plays the role of an effective soft parti¬ 
cle diameter. For distances r > a s , two Hertz model par¬ 
ticles do not interact with each other. The strength of the 
continuous potential is quantified by the non-dimensional 
elasticity parameter (effective potential strength), 


2 ^ 

15k B T(l-is 2 ) ’ 


( 2 ) 


depending on the bulk modulus Y and the Poisson ratio 
v of a particle [53, [H|. On approximating the mesoscopic 
particle elastic moduli Y and v by macroscopic values, 
the estimate e « 10 4 ~ 10 5 is obtained for micron-sized 
microgels. Note here the strong size dependence e oc of, 
implying a significantly decreased potential strength for 
smaller microgel particles. The Hertz potential has been 
shown to provide a nearly fit-parameter free description 
of the structure and phase behavior of neutral micro¬ 
gel systems, in good agreement with experimental results 
for various particle sizes [53, EH , and values of e in the 
range from 10 2 — 10 4 . While for large e the Hertz and 
hard-sphere potentials lead to very similar results for the 
equilibrium microstructure, the effective diameter in the 
Hertz potential is in general somewhat larger than the 
corresponding hard-sphere value, i.e. a s > a. This re¬ 
flects the fact that the Hertz potential incorporates over¬ 
all the softness of a microgel particle which in turn orig¬ 
inates from the radially inhomogeneous crosslinker den¬ 
sity. 

For very large values of e, the Hertz potential is prac¬ 
tically indistinguishable from the hard-sphere potential 


rr / \ Joo r < a = 2a 
V m (r) = < (3) 

10 r > a , 

with o s playing now the role of the hard-sphere diameter 
a. Considering the high-energy deformation penalty e of 
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micron-sized, non-ionic PNiPAM microgels, the usage of 
the hard-sphere potential in place of the Hertz potential 
is well justified (see inset in Fig. [3]). This is advantageous 
from a theoretical viewpoint since colloidal hard spheres 
are among the most thoroughly studied soft matter sys¬ 
tems. 

In addition to being the key quantities in static scatter¬ 
ing experiments on concentrated suspensions, the static 
structure factor, S(q), and its associated RDF, g(r), 
are required as inputs in various methods of calculating 
colloidal transport properties. For the Hertz potential 
model, we determine the two pair correlation functions 
numerically from solvin g the ap proximate Percus-Yevick 
(PY) integral equation [45], [46| . The corresponding pair 
correlation functions for the hard-sphere model are de¬ 
termined from the analytic PY solution combined with 
the Verlet-Weis (VW) correction [47| incorporating the 
accurate Carnahan-Starling equation of state (see also 
HI)- The VW correction compensates in particular the 
overestimation by the PY solution of the principal peak 
height, S(q m ), of the hard-sphere structure factor for vol¬ 
ume fractions (f> = ( ir/6)na 3 > 0.4. Here, n is the number 
concentration of particles. 



FIG. 1. Comparison of the experimental static structure fac¬ 
tor, S(q), of PNiPAM microgels (filled symbols, taken from 
[i3|) with the Verlet-Weis corrected Percus Yevick prediction 
(solid lines), for various particle volume fractions 4> as indi¬ 
cated. The bare PY structure factor for the largest considered 
volume fraction (f> = 0.5 is represented by the dashed curve. 

In Fig. [T] we compare the PY-VW structre factors for 
the hard-sphere model with static light scattering results 
by Eckert and Richtering [43| for a concentration series 
of non-ionic poly(N-isopropylacrylamide) (PNiPAM) mi¬ 
crogels in DMF. The details of the particle synthesis are 
given in [H, [49|, [50|- In our PY-VW calculations, we 
have used the experimentally obtained (mean) particle 
diameter cr = 240 nm, and the (volume-swelling cor¬ 
rected) volume fractions cf> as given in [43]. We refer 
to this reference for the details on how the volume frac¬ 


tions have been determined. The agreement between the 
experimental and PY-VW structure factors is good even 
for small wavenumbers q. There is in particular a sig¬ 
nificantly improved agreement for the largest considered 
volume fraction (dashed curve in Fig. [I]), as compared 
to the bare PY S(q) used in the earlier work [43|. The 
remaining small deviations for the most concentrated sys¬ 
tem at (j) = 0.5 can be at least partially attributed to ex¬ 
perimental uncertainties, and possibly also to the break¬ 
down of the exact factorization of the mean scattered 
intensity into static structure and form factors as dis¬ 
cussed by Likos et al. [5l|, [52}. In Fig. (2J the experimental 



FIG. 2. Experimental structure factor at 4> = 0.378 (filled 
circles) in comparison with the best-fit hard-sphere model PY- 
VW S(q) (black solid line), and the best-fit Hertz model S(q) 
calculated in PY approximation (red solid line). For the Hertz 
model, the parameters e = 10 4 , <r a = 245 nm, and (p s = 0.398 
have been used. 

structure factor of PNiPAM microgels is compared, for 
(j) = 0.378, with the best-fit hard-sphere and the best-fit 
Hertz potential structure factors. For the Hertz model 
system, the parameter values a s = 245 nm, e = 10 4 , and 
(j> s = (7t/ 6)nof = 0.398 have been used. The theoretical 
curves coincide practically, as expected for the invoked 
large value of e which restricts the softness range of the 
Hertz potential to a narrow interval around the effec¬ 
tive diameter a s . This can be noticed from the inset of 
Fig. [3] The remnant potential softness still necessitates a 
slightly larger effective diameter cr s than the hard-sphere 
model value a = 240 nm, for the Hertz potential incor¬ 
porates the interactions by dangling polymer chains at 
the periphery of microgel particles. The larger particle 
size in the Hertz potential model goes along with a larger 
volume fraction cj) s than that used in the best-fit hard- 
sphere model result. With increasing volume fraction 
in the range </> s £ [0.299 — 0.5], a slight decrease of the 
best-fit effective diameter cr s is observed. However, ow¬ 
ing to the smallness of the particle size shrinkage, for the 
considered PNiPAM in DMF system this effect can be 
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disregarded in its influence on static and dynamic prop¬ 
erties. 

The g(r)’s corresponding to the best-fit Hertz and hard- 



FIG. 3. Radial distribution functions, g(r), corresponding 
to the hard-sphere and Hertz potential structure factors in 
Fig. [2] The inset compares the hard-sphere and Hertz pair 
potential, for the parameters in Fig. [2] 

sphere model structure factors in Fig. [2] are depicted 
in Fig. [3] The slight softness of the Hertz potential is 
reflected in the somewhat reduced peak height of g(r) 
(red solid curve), and in the slight extension of the RDF 
into the overlap region r < a s . From Figs. [2] and [31 
we conclude that the pair correlations of the considered 
PNiPAM suspensions in the experimental concentration 
range are fully compatible with the simple hard-sphere 
potential. For microgels of lower crosslinker density and 
smaller size where the softness of the particles is impor¬ 
tant, the Hertz potential model should be used. 


B. Hydrodynamic Particle Modeling 

As an illustration of the concept of a hydrodynamic 
radius ah and its associated slip length Lh , consider two 
simple hydrodynamic particle models, namely uniformly 
fluid-permeable rigid spheres [H,[54|], and non-permeable 
rigid spheres with Navier par tial slip hydrodynamic sur¬ 
face boundary condition [55| . The physical (material) 
sphere radius in both models is denoted here by a. 

In the uniformly-permeable sphere model, the pore- 
size averaged fluid flow inside a particle is described by 
the Brinkman-Debye-Bueche (BDB) equation, and the 
outside flow by the low-Reynolds-number Stokes equa¬ 
tion [5f|. The imposed boundary condition is here that 
the fluid velocity and tangential stress are continuous 
across the particle surface. The suspension transport 
properties in this model depend on the material-specific 
parameter 



equal to the ratio of the hydrodynamic penetration 
length, 1 /k, and particle radius a. The penetration 
length is roughly equal to the mean pore size of the rigid 
skeleton of a permeable sphere. In the limit X x —> 0 of 
vanishing mean pore size, a non-permeable sphere with 
no-slip hydrodynamic boundary condition (BC) on its 
surface is obtained. For the BDB equation to apply, the 
mean pore size should be no larger than one tenth of the 
particle radius so that X x < 0.1. The present model has 
been generalized to spherical particles with a permeabil- 
ity profile n(r) varying with the radial distance (see, e.g. 
[5ij |) such as core-shell particles [2ll l57l - l59| consisting of 
a dry core and a permeable outer shell. The core-shell 
model describes in a coarse-grained manner the hydro- 
dynamic effect, e.g., of a polymer brush surrounding the 
core. 

The Navier partial-slip model on the other hand de¬ 
scribes fluid-impermeable colloidal spheres where the 
fluid is allowed to partially slip along their surfaces. The 
associated Navier partial-slip BC demands for a station¬ 
ary sphere at each surface point the proportionality of 
surface-tangential fluid velocity, un, and shear-stress, t m , 
according to 



( 5 ) 


The proportionality constant is given by the ratio of the 
so-called Navier length, l s , and the fluid shear viscos¬ 
ity rjQ. In the limit l* —>• 0, the no-slip BC describing 
zero surface slip is recovered. Here, l* = l s /a is the re¬ 
duced Navier length. In the opposite limit l* —» cso, the 
free-surface BC of zero tangential stress is obtained, cor¬ 
responding to fluid perfectly slipping along the sphere 
surface in form of local plug flow. The Navier partial- 
slip BC can serve as an effective description of a hy¬ 
drophobic particle surface, and of a rigid particle with 
surface roughness and corrugations (60|. It is also appli¬ 
cable when non-adsorbing (short) polymers are dispersed 
in the fluid, owing to the formation of a thin clear-fluid 
depletion layer at the particle surfaces [ 6 ll ]. 

The hydrodynamic radius is a single-particle property 
which depends on the intra-particle hydrodynamic struc¬ 
ture of the spherical particle, and in principle also on the 
considered single-particle transport property. It can be 
defined operationally through the Stokes-Einstein-Debye 
expressions 


Dl(a) 

D r 0 (a) 


k B T 

a^(a) 

k B T 

8mi 0 a r h (a) 3 ’ 


( 6 ) 

( 7 ) 


for the translational and rotational diffusion coefficients 
Dq and Dq, respectively, of an isolated spherical particle 
in an infinite fluid. An additional definition not consid¬ 
ered here is based on the intrinsic viscosity, [ 77 ], of dis¬ 
persed spheres 0 . Here, a denotes a set of parameters 
characterizing the hydrodynamic particle red structure. 
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For the two considered models is a = {X x ,l*}. The re- 
spective translational and rotational hydrodynamic radii 
are (5/ 


a h(^x) 

a 

a 


2x 2 (x — tanh (a;)) 

2a ; 3 + 3(a; — tanh (a;)) 

3 3 coth (a;) 1 3 


( 8 ) 

(9) 


for a homogeneously permeable sphere with x = 1/X x , 
and [55j 


a 

a r h(lt) 

a 


1 + 2 1 * 

1 + 3 1 * 
i 

1 + 31*) 


( 10 ) 

( 11 ) 


for a Navier partial-slip sphere. 

One can associate each hydrodynamic radius ah with 
a reduced hydrodynamic slip length, L\, through 


L* h (a) 


Lh{a) _ 1 _ ah(u) 
a a 


( 12 ) 


The quantity L ^ is the relative radial distance, Lh = 
a — ah, of the apparent no-slip spherical surface inside 
the particle to its outside surface, in units of the out¬ 
side surface radius a. Except for X x , the asterisk is used 
throughout to label lengths given in units of a, . 

In the experimentally common situation where L £ is 
small, curvature effects are negligible and the particle- 
fluid interface can be described as flat. As it is explained 
in [23j |. the key point to notice is that the reduced slip 
length, L* h j = 1 ~ ahj/a , in the flat-interface approx¬ 
imation and its associated flat-interface hydrodynamic 
radius, ahj , are independent of the single-particle trans¬ 
port properties D q, Dq and [ rj ] used in their definition. 
Since each of these transport properties is associated with 
a particular ambient velocity field, e.g. with a constant 
flow field in case of Ug, an equivalent statement is that 
Lhj is independent of the ambient flow. 

It was shown in [23| that the relations 

L* h (a)=L* hJ (a)[l + 0(L* hJ (a))] (13) 

a* h (a)=a* hJ (a)[l + 0(L* hJ (a) 2 )] , (14) 

apply universally to rigid spheres with arbitrary, non¬ 
singular spherically symmetric hydrodynamic structures 
and boundary conditions, independent of the invoked 
single-particle transport property. 

These relations are readily verified, using Eqs. © - 
CD, for the special case of a Navier partial-slip sphere 
where the flat-interface slip length is equal to 

L*hj(l* s )=l*s, ( 15 ) 

and for a uniformly permeable sphere where 


Note that the slip length Lh is equal to the Navier length 
l s in the flat-interface limit only. 

Eqs. (11311141) do not apply to non-rigid particles such 
as spherical liquid droplets, and rigid particles of singular 
hydrodynamic structure. The latter case is exemplified in 
|5J] by a hollow sphere with a uniformly permeable rigid 
shell of infinitesimal thickness. The hydrodynamic radii 
a\ and a r h of the ultra-thin hollow sphere differ already 
to linear order in the smallness parameter A^. As for a 
droplet, an inscribed apparent no-slip spherical surface 
which is necessarily rigid can not be introduced for the 
singular hollow-sphere example. 

Eqs. (l8lfTTl) can be inverted to obtain the material- 
specific parameters X x and l* in terms of the ratio 


or likewise in terms of the reduced width (slip length), 
7 = 1 — 7 , of the apparent fluid annulus region sur¬ 
rounding the apparent no-slip sphere. The parame¬ 
ter 7 completely characterizes the intra-particle hydro- 
dynamic structure in the hydrodynamic radius model 
(HRM) where spherical particles of arbitrary hydrody¬ 
namic structure are described hydrodynamically as no- 
slip spheres of reduced radius a/j, for unchanged pair po¬ 
tential. This implies in particular an unchanged excluded 
volume radius a > ah- 

Numerical results illustrating this inversion are de¬ 
picted in Fig. U where the reduced penetration length of 
a permeable sphere, A^’ r , and the reduced slip length of a 
Navier partial-slip sphere, (l*) t,r , are plotted as functions 
of 7 . The reduced lengths derived from the translational 
(rotational) Stokes-Einstein- (Debye) relation are labeled 
by the superscript t (superscript r). For a thin annulus 



FIG. 4. Reduced fluid penetration length X x r (red lines) and 
reduced Navier length (l*) t,r (blue lines), as functions of the 
reduced fluid annulus width 7=1 — 7 . Solid lines (dashed 
lines) describe quantities derived from the translational (ro¬ 
tational) single-sphere Stokes-Einstein-(Debye) relation. 


Lhj(^x) ~ 


(16) shell with 7 > 0.9, the four reduced lengths are small and 
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commonly described by the dashed-dotted line 7 , corre¬ 
sponding to \ x « Lh and l* ss Lh- This situation is 
common to colloidal systems for which values < 0.05 
are typically found, and for which the Navier length is 
commonly equal to a few nanometers. 

Differences between the considered lengths are signif¬ 
icant for larger values of 7 . For a given reduced hydro- 
dynamic radius 7 = ah/a, two different values Aj, ^ \ r x 
are then obtained. The perfect-slip limit l* —>• 00 cor¬ 
responds to 7 = 1/3 for translational, and to 7 = 1 for 
rotational diffusion. While a rotating perfect-slip sphere 
experiences no friction, akin to a point particle for which 
7 = 0, the values of 7 for a translating Navier sphere 
do not exceed 1/3. Regarding the translational diffusion 
coefficient, a perfect-slip sphere with 7 = 2/3 has the 
same value, Dq = kBT/(4nr]oa), as a uniformly perme¬ 
able sphere with A^, ss 0.278. However, such a large value 
of \ x is unrealistic. In order for the BDB equation to de¬ 
scribe flow inside a permeable particle, X x should be no 
larger than 0 . 1 , corresponding to a maximally allowed 
penetration length equal to one-tenth of the particle ra¬ 
dius. 


III. SHORT-TIME TRANSPORT PROPERTIES 

In studying the dynamics in dispersions of colloidal 
particles, one distinguishes the short-time regime where 
tb <C t <C T£> from the long-time regime where t tb- 
Here, t denotes the correlation time of particle concen¬ 
tration fluctuations, and tb and td are the character¬ 
istic particle momentum relaxation and diffusion times, 
respectively. The time tb is the decay time of velocity 
correlations of a particle. For submicron sized particles 
in a low-viscosity fluid such as the here considered PNi- 
PAM in DMF microgels, tb is of the order of microsec¬ 
onds. The particle diffusion time can be estimated by 
tb ~ a 2 /Dq which for the considered microgels is in the 
millisecond range. In the short-time regime, the particle 
configuration has changed so little that the slowing influ¬ 
ence of non-hydrodynamic direct interactions is not yet 
operative, different from the solvent-mediated His which 
act quasi-instantaneously on the colloidal scale. Short- 
time transport properties are thus expressible as simple 
equilibrium averages, with the direct interactions enter¬ 
ing only by their influence on the equilibrium microstruc¬ 
ture, e.g on the RDF g(r). Long-time transport proper¬ 
ties are additionally influenced by the Brownian motion 
of the particles so that they depend in a direct way both 
on direct and hydrodynamic interactions. 

In this section, we describe our analytic toolbox of cal¬ 
culating short-time diffusion and rheological properties of 
hydrodynamically structured particles, in the framework 
of the HRM. In the subsequent sections, the toolbox is 
used for analyzing light scattering data on PNiPAM in 
DMF suspensions. 


A. Hydrodynamic function 

The colloidal short-time regime is probed in DLS ex¬ 
periments by analyzing the initially exponential decay, 

S(q,t <C t d ) = S/g/exp [~q 2 D(q)t] , (18) 

of the dynamic structure factor, S(q,t), in its depen¬ 
dence on the modulus q of the single-scattering wave vec¬ 
tor q. The dynamic structure factor is the g-th Fourier 
component of spatio-temporal correlations in thermally 
induced particle concentration fluctuations. On basis 
of the many-particlc generalized Smoluchowski equation 
(GSmE), the wavenumber-dependent diffusion function, 
D(q ), quantifying the short-time exponential decay can 
be expressed by the ratio [HJ 
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of the hydrodynamic function H ( q ) and the static struc¬ 
ture factor S(q). The hydrodynamic function is the key 
quantity containing information on colloidal short-time 
diffusion processes. It can be expressed by the equilib¬ 
rium configurational average 0, 

= Eq-^ (x) -qe iq[Ri - R ^) , (20) 

over an ensemble of N colloidal spheres, with the thermo¬ 
dynamic limit N — > 00 for a fixed particle concentration 
n performed subsequently, in order to describe a macro- 
scopically large scattering volume. Here, q is the unit 
vector in the direction of q, and ksT is the system tem¬ 
perature T multiplied by the Boltzmann constant. The 
translational mobility tensor /x y - linearly relates the hy¬ 
drodynamic force acting on particle j at the instant posi¬ 
tion Rj to the resulting velocity change of a particle i. It 
depends on the hydrodynamic structure of the globular 
particles, e.g. on particle sizes, surface BCs and fluid per¬ 
meability, and on the many-particles His for an instant 
configuration, X = {Ri, • • • , Rat}, of particle centers. 

Without His, H(q) would be identically equal to one 
independent of q. Any wavenumber dependence reflects 
thus the presence of His. The hydrodynamic function 
consists of a g-independent self-part and a distinct part, 
i.e. 

H(q) = ^ + H d (q), (21) 

with the distinct part, H d (q), approaching its asymptotic 
value zero for q —> 00 . Thus, for large q corresponding to 
small distances resolved, H(q) reduces to the short-time 
self-diffusion coefficient, Ds, expressed in units of the 
single-sphere translational diffusion coefficient D q. The 
coefficient Ds quantifies the initial (i.e., short-time) slope 
of the mean-squared displacement of a particle in pres¬ 
ence of other ones. Owing to the slowing influence of the 




His at non-zero particle concentrations, Dg is smaller 
than Dq to which it reduces at infinite dilution only. 

The hydrodynamic function has the physical interpre¬ 
tation of a short-time generalized sedimentation coeffi¬ 
cient for a homogeneous suspension of monodisperse col¬ 
loidal spheres subjected to a weak force field colinear with 
q and oscillating spatially as cos (q • r) with positions r 
[62 ]. For q —> 0, a homogeneous force field is recovered 
such as the local gravitational field on earth. Conse¬ 
quently, 


K=-^ = ]imH(q) (22) 

Vo ?->o 

is the short-time average sedimentation velocity, I4ed, of 
hydrodynamically interacting colloidal particles, normal¬ 
ized by the particle model dependent mean velocity, Vo, 
of an isolated particle sedimenting in the same force field. 
The sedimentation coefficient K is related to the short- 
time collective diffusion coefficient, Dc, by 

= (23) 

where for monodisperse particles 5(0) = Hm q ^oS(q) is 
the relative osmotic compressibility of the suspension. 
At larger concentrations where 5(0) is small, Dc is sig¬ 
nificantly larger than Dq. In principle, the short-time 
coefficient Dc should be distinguished from the associ¬ 
ated long-time collective diffusion coefficient appearing 
in Fick’s law of macroscopic gradient diffusion. The lat¬ 
ter, however, is only slightly smaller than the short-time 
coefficient, by less than 7% even for a concentrated sus¬ 
pension of no-slip colloidal hard spheres at volume frac¬ 
tion (f> = 0.45 [63|. The difference between the two col¬ 
lective diffusion coefficients can be expected to be even 
smaller for particles with weaker His such as perme¬ 
able and partial-slip spheres, and charge-stabilized par¬ 
ticles with long-range electrostatic repulsion. Different 
from collective diffusion and sedimentation, the long-time 
translational self-diffusion coefficient Dl can be substan¬ 
tially smaller than its short-time sibling Dg, owing to 
the retarding relaxation of next-neighbor particle cages 
which are slightly distorted from their equilibrium spheri¬ 
cal symmetry at long times. For colloidal hard spheres at 
the freezing transition concentration, e.g., Dg/Dg « 0.1 
accord ing to the empirical freezing rule by Lowen et 

alMm 

In this work, we discuss useful analytic scaling rela¬ 
tions for the calculation of short-time diffusion and vis¬ 
cosity properties of suspensions of internally structured 
particles. These relations are complemented by two semi- 
analytic and easy-to-implement approximation schemes 
which have been successfully applied in the past to sus¬ 
pensions of no-slip neutral and charge-stabilized spher¬ 
ical particles [42], I66l - l68| . The first scheme is the so- 
called pairwise-additivity (PA) approximation where all 
two-body contributions to the hydrodynamic mobilities 
fiij are accounted for including near-contact lubrication 


terms, but three-body and higher-order contributions are 
disregarded. The PA method is thus well suited for low 
concentrations where it is basically exact. The second 
method is a self-part corrected improved version [42j of 
the Sj renormalized concentration fluctuation expansion 
by Beenakker and Mazur [69j - l7l| . where many-particle 
His contributions are approximately accounted for in 
terms of so-called ring diagrams, but with lubrication 
corrections disregarded (see also [72]] ) ■ The dy method is 
well suited in particular for larger concentrations. The 
only input required by these methods is the RDF or like¬ 
wise the static structure factor. The explicit expressions 
by both methods for H (q) and the high-frequency (short- 
time) suspension viscosity r )(*>, generalized to the HRM, 
are presented in the Appendices [A] and [B] respectively. 

A key observation regarding the translational mobil¬ 
ity tensors fiij of hydrodynamically structured spherical 
particles is the relation 




1 + 0(L* h /) 


(24) 


where /l^hrm (X; a-hj) are the mobility tensors of the 
associated HRM. The relation follows from Eqs. m 
and (m in conjunction with a general scattering series 
expansion of the exact IV-sphere translational mobility 
tensors [22], [23| . Since the short-time transport proper¬ 
ties are equilibrium averages of specific mobility tensor 
elements, it follows that 


H(q) = Hhrm(?) 


1 + 0(L* h f) 


(25) 


with an analogous expression valid for r^oo. Thus, the er¬ 
ror introduced in calculating short-time transport prop¬ 
erties of hydrodynamically structured particles using the 
simplifying HRM is of 0(L^ j 2 ) small. 


B. Translational self-diffusion 

In [341 Iril , it was shown by comparison with computer 
simulation results for uniformly permeable hard spheres 
with X x < 0.2 that Dg normalized by its infinite dilution 
Dq is well represented by the simple scaling relation, 


^TPTTr = 1 + At <£ (1 + 0-120 - 0.70(/) 2 ) , (26) 

L>Q (A t J 

for all </> < 0.5 and with an accuracy of 3.5% or better. A 
similar scaling relation has been obtained for the short- 
time rotational self-diffusion coefficient of permeable- 
sphere suspensions with hard-core interactions [HI] . The 
only dependence in Eq. (ESI) on the permeability param¬ 
eter A* characterizing the intra-particle hydrodynamic 
structure is contained in the second virial coefficient, 
A t = X t (\ x ), which is the linear coefficient in the expan¬ 
sion of Dg in powers of cj). Eq. (ESI) states that the Dg 
for hydrodynamically structured particles can be scaled 
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to the corresponding coefficient of no-slip hard spheres 
where = 0 and = 0) = —1.8315. 

In fact, on the level of the HRM model Eq. (l26l) is ap¬ 
plicable to suspensions of hard spheres of arbitrary hy¬ 
drodynamic structure, with deviations from the Ds of 
the actual particle system being of quadratic order small 
in the reduced slip length L £ according to Eq. C5l) . On 
recalling that for spherical particles with hard-core di¬ 
rect interactions only the HRM reduces to the spherical 
annulus model, the only input required in Eq. (12611 is 
the second virial coefficient, At ( 7 ), of spherical annulus 
particles depending on the ratio, 7 = ah/a, of hydro- 
dynamic and hard-core radius. Using for At ( 7 ) the Eq. 
(IA 6 I) in Appendix [A] in conjunction with tabulated nu¬ 
merical values for the longitudinal and transversal mobil¬ 
ity coefficients, in(r) and yn(r), of two no-slip spheres 
of radius ah calculated using the method of Jeffrey and 
Onishi [73| . we have obtained numerically precise values 
for the second virial coefficient of spherical annulus par¬ 
ticles. These values are described to high accuracy by 
the polynomial fit 

A t ( 7 ) = - 1.8315 + 7.8207 - 14 . 2317 2 (27) 

+ 14.9087 3 - 9.383T 4 + 2.7177 s , 

accounting for the numerically correct limiting value 
A t ( 7 = 1 ) = —1.8315 of no-slip hard spheres. 



FIG. 5. Modulus |At | of the first-order virial coefficient of Ds 
in the spherical annulus model, as a function of the reduced 
fluid annulus shell width 7 = 1 — 7 . Solid line: Polynomial 
fit according to Eq. (1271) . Filled squares: Tabulated values 
by Cichocki et al. [22|] for spherical annulus particles with 
7 < 1/3. Filled diamonds: Tabulated values by Cichocki and 
Felderhof ItI . 

In Fig. [5j the high accuracy of the polynomial fit in 
Eq. (1271) is shown in comparison with earlier numerical 
results [ 22 I. |74] for the second virial coefficient of spheri¬ 
cal annulus particles. Note that A t (7 = 0) = 0 relates to 
the limiting case of spherical annulus particles interact¬ 
ing hydrodynamically as point particles for which there 


are no hydrodynamic self-reflections so that Ds reduces 
to Dq independent of (j>. The associated long-time coeffi¬ 
cient Dl for 7 = 0 is still ^-dependent and smaller than 

n* 

u o- 



FIG. 6 . Results for the reduced self-diffusion coefficient, 
Ds/Dq, of the spherical annulus model as function of 7 = 
1 — 7 , for several volume fractions 4> as indicated. Solid lines: 
Prediction by the scaling formula in Eq. © in combination 
with the first order virial coefficient fitting polynomial in Eq. 
(1271) . Closed symbols: HYDROMULTIPOLE simulation data 
for the spherical annulus model tabulated in [21] . Open sym¬ 
bols: Simulation data for uniformly permeable spheres tabu¬ 
lated in na. with Eq. © used in mapping the permeability 
parameter A x onto the annulus model parameter 7 = ah/a 
(recall Fig. ©. 

The accuracy of the analytic scaling formula in Eq. 
(ESI) in combination with Eq. 03 for the short-time 
self-diffusion coefficient of spherical annulus particle sys¬ 
tems is established in Fig. [ 6 ] by the comparison with 
high-precision benchmark simulation data for the spheri¬ 
cal annulus [Til and uniformly permeable particle models 
E3- Regarding the latter model, the conversion of \ x to 
the related reduced hydrodynamic radius parameter 7 in 
the spherical annulus model was done using Eq. J 8 |) for 
Dq(A x ). This corresponds to the inversion of the curve 
for A* in Fig. (©) in terms of 7 . For example, the small¬ 
est reasonably selected value A^ = 0.1 corresponds to 
7 = 0.89. Results for Ds are depicted in Fig. [ 6 ]in depen¬ 
dence on the reduced slip length 7 , for volume fractions 
extending over a broad concentration range. 

The excellent agreement between the scaling formula 
for Ds (solid lines) and the simulation data (symbols) 
does not only validate this formula. For the special case 
of permeable hard spheres, it additionally highlights the 
good performance of the HRM, for a broad concentration 
range and permeability values largely exceeding the ones 
discussed earlier in the thin shell-limit discussion of the 
core-shell model in [ 22 j. 

In summary, the analytic formula in Eqs. (1261) and (1271) 
allows for a quick and accurate calculation of the trans- 
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lational short-time self-diffusion coefficient of hydrody- 
namically structured spherical particles with hard-core 
interactions. The only input is the single-particle prop¬ 
erty ah/a which can be determined experimentally, e.g., 
by a DLS measurement of D q in conjunction with a static 
scattering experiment determining the excluded volume 
radius. The formula for Dg is applicable also to particles 
with a short-range, weakly soft pair potential such as the 
Hertz potential for non-small potential strengths e, pro¬ 
vided the effective diameter a s and the related volume 
fraction <f> s are used instead of a and <fi. 

The formula for Dg in Eqs. (ESI) and m based on 
its virial expansion is not valid for particles with a long- 
range, soft pair potential where an (effective) excluded 
volume diameter is not the characteristic parameter. An 
example in case are low-salinity suspensions of charge- 
stabilized colloidal spheres, where according to theory, 
simulation and experiment Dg has to good accuracy 
the fractional concentration dependence Dg/D° s — 1 « 
Ag 0 4 / 3 , with a coefficient Ag ~ 2.5 — 2.9 which varies to 
a small extent with the particle size and charge [67, 75} . 
The initial slope of Dg at 0 = 0 is thus zero in these 
charge-stabilized systems. Different from the formula in 
Eqs. (l26l) and (|27l) . the HRM is applicable also to spher¬ 
ical particles with arbitrary soft direct interactions, and 
to particles of fuzzy hydrodynamic structure without a 
sharp outer boundary. In the framework of the HRM, 
the coefficient Dg and other short-time transport prop¬ 
erties of particles with soft interactions can be approx¬ 
imately and semi-analytically calculated using the PA 
method at smaller and the self-part corrected Ay method 
at larger concentrations. The inputs g(r) and S(q) to 
these methods can be obtained from Ornstein-Zernike 
integral equation schemes such as the analytical scheme 
for charge-stabilized particles introduced in [zE Q. 


C. Sedimentation coefficient 


Different from self-diffusion, the neglect of hydrody¬ 
namic near-distance and particularly lubrication effects 
in approximate analytic calculations is less consequential 
regarding sedimentation. Therefore, A can be described 
semi-qualitatively on the Rotne-Prager (RP) type level 
where only the long-distance dipolar contribution to the 
mobility tensors is accounted for. This amounts hydro- 
dynamically to the neglect of all hydrodynamic flow re¬ 
flection by the particles. 

Using the RP approximation in conjunction with the 
analytic PY solution for the Laplace transform of rg(r), 
Contreras-Aburto et al. [78} have derived analytic ex¬ 
pressions for the short-time sedimentation coefficient of 
Navier partial-slip and uniformly permeable spheres with 
hard-core interaction. Here, we present the according ex¬ 
pression for spherical annulus particles, 


A'rp (4>, 7 ) 


1 +')<j) 



H 2 - 0) - io ]\ 

20(l + 2<£) \) 


, (28) 



FIG. 7. Concentration dependence of the sedimentation co¬ 
efficient, A, of the spherical annulus model for reduced hy¬ 
drodynamic radius values 7 as indicated. Solid lines: RP 
approximation Arp in Eq. (1281) . Filled symbols: Simulation 
data taken from [ 2 lj . 


which includes hydrodynamic point particles as a limit¬ 
ing case for which A(0, 7 = 0) = 1. In the opposite limit, 
7 = 1, of no-slip hard spheres, E q. (E51) reproduces an ex¬ 
pression by Banchio and Nagele [67] which was rederived 
subsequently by Gilleland et al. [79( using a variational 
method. 

The comparison in Fig. [7] of the RP based analytic 
formula in Eq. (1281) with benchmark simulation data for 
spherical annulus hard spheres J 2 l[ shows that the sed¬ 
imentation velocity is overestimated by the formula at 
larger concentrations. This is a consequence of the ne¬ 
glect of flow reflections in the RP approximation which 
becomes less severe with decreasing annulus parameter 7 , 
owing to the for a fixed 0 increasing distances between 
the hydrodynamic particle surfaces. For the lowest con¬ 
sidered value 7 = 0 . 6 , excellent agreement between the 
simulation data and Arp is observed. The largest devia¬ 
tions occur for no-slip hard spheres where Arp provides 
an upper bound to the exact sedimentation coefficient 
[62, [79j. As an aside, we note that even at 0 = 0.5, Arp 
changes only slightly if the VW-corrected g(r) is used 
instead of the bare PY g(r). 

While Arp nicely describes the trends of the exact 
spherical annulus sedimentation coefficient A in its cf> and 
7 dependence, in search of an improved analytic expres¬ 
sion we make the ansatz 

A ( 0 , 7 ) = 1 + A k ( 7 ) Mk ( 0 , 7 ) ( 29 ) 

= 1 + Ak ( 7 ) 0 [1 + O (</>)] , 

where Ak ( 7 ) is the first-order virial coefficient of the sedi¬ 
mentation coefficient for which high-precision values have 
been provided by Cichocki et al. [22, M- According to 
Fig. [51 these tabulated values are well represented by the 
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FIG. 8 . First-order virial coefficient, A k, of the sedimenta¬ 
tion coefficient of spherical annulus particles. Filled symbols: 
Tabulated values by Cichocki et al. for thick (diamonds) [zH 
and thin annulus shell (squares) systems ( 22 j, in comparison 
with the polynomial in Eq. (1301) (black solid line). 


FIG. 9. Concentration dependence of the scaling function 
mk (< i >, 7 ) associated with the sedimentation coefficient of 
spherical annulus particles, for values of 7 as indicated. Col¬ 
ored closed symbols: Values obtained from simulation data of 
If ( 0 , 7 ) [ 21 ]. Colored solid lines: Semi-empirical formula in 
Eq. (13T1) . 


polynomial, 

(7) = 


-6.5464 + 8.5927 - 3.9017 2 
+ 2 .OII 7 3 - 0 . 1427 4 


(30) 


in the full parameter range 0 < 7 < 1. At 7 = 1 in 
particular, the known value Ak = —6.5464 of no-slip hard 
spheres is recovered from the polynomial. 

In the first equality in Eq. (l29l) . we have introduced 
the sedimentation scaling function uk = (K — 1 ) /Ak- In 
34J, a scaling ansatz analogous to Eq. (ITil) was used for 
the short-time translational and rotational self-diffusion 
coefficients of uniformly permeable hard spheres. By the 
comparison with simulation data, it was shown therein 
that the scaling functions us{4>, X x ) and ur( 0, A x ) asso¬ 
ciated with translational and rotational self-diffusion, re¬ 
spectively, are practically independent of the permeabil¬ 
ity coefficient for all values X x < 0.1. They are therefore 
well approximated by (i.e. scaled to) the respective func¬ 
tions ms,r(0, X x = 0) of non-permeable solid spheres. On 
using a third-order polynomial fit of us(4>,X x = 0) ob¬ 
tained from simulation data of no-slip hard spheres, and 
the first two known virial coefficients of X x = 0), 

Eq. (l26l) for Ds(4>,X x ) has been obtained [34 142] . 

As noticed already in the context of permeable spheres 
a uk depends significantly on the intra-particle hydro- 
dynamic structure, different from its self-diffusion sib¬ 
lings. In Fig. H this is demonstrated for the spheri¬ 
cal annulus system using tabulated simulation data for 
7 ) in [2ll| . The simulation data of uk(^ 7 ) for a 
fixed (f> > 0 clearly differ from each other for different 
7 values. Thus, the 7 -dependence of K(cf>, 7 ) cannot be 
embodied solely in terms of the first-order virial coeffi¬ 
cient, in contrast to Eq. (1261) describing Ds . However, 
as it is shown in Fig. [9] the simulation data for uk are 


well described for 7 < 0.8 by a forth-order polynomial in 
7 (f>, namely 


UK (<(>,7) = 01 


3.3487 <t> + 7.426 ( 7 (j)f 


-10.034 (rfc/)) 3 + 5.882 ( 7 0) 4 


(31) 


The resulting analytic expression, 

If(0,7) = l- A K (7)0[1-3.34870 (32) 

+7.426 ( 7^) 2 - 10.034 ( 7 </>) 3 + 5.882 ( 7 0) 4 ] , 


for the sedimentation coefficient provides in conjunction 
with Eq. (I3U1) for Xxi'y) an accurate description in the 
from the experimental viewpoint sufficiently broad pa¬ 
rameter range 7 £ (0.8 — 1}. The numerical coefficient 
3.348 in the bracket of Eq. (El is selected such that at 
7 = 1 the correct numerical value 21.918 of the second 
virial coefficient of no-slip rigid spheres [74[ is recovered. 


From Fig. [TO] the good agreement of the semi- 
empirical formula for K(<f>, 7 ) in Eq. (ITil) with the spher¬ 
ical annulus simulation data [|I| is noticed for 7 > 0 . 8 . 
The figure depicts also simulation data for permeable 
spheres where X x has been converted to the respective 
7 using Eq. ©■ Moreover, results for 7 ) are dis¬ 
played as predicted by the self-part corrected J 7 scheme 
with VW-PY structure factor input where the self-part 
contribution to K was calculated according to Eq. El- 
See Appendix iBl for details on the £7 method of calcu¬ 
lating H(q) and its self-part correction. While in good 
overall agreement with the simulation data, the self-part 
corrected £7 method results for K deviate significantly 
at larger (j>. Different from Eq. (13^1) which applies to 
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FIG. 10. Sedimentation coefficient of the spherical annulus 
model as a function of <j), for values of 7 as indicated. Filled 
squares: Simulation data for spherical annulus particles [2l| . 
Solid lines: Analytic formula in Eq. (1321) with Ak according 
to Eq. (13011 . Dashed lines: Self-part corrected S'y method 
prediction, with self-part Ds according to Eq. (HD and VW- 
PY input for S(q). 


FIG. 11. Reduced hydrodynamic function, hd(q), for spherical 
annulus particles at fixed volume fraction 0 = 0.35 and vary¬ 
ing 7 as indicated. The wavenumber is scaled in terms of the 
hard-sphere diameter a = 2a. Solid lines: 5 7 method results 
using VW-PY structure factors as input (cf. Appendix [B}. 
Filled squares: Simulation result for no-slip (ns) hard spheres 
( 7 = 1 ) taken from flsil . 


spherical particles with hard-core interactions only, the 
5"/ method is applicable also to particles with soft inter¬ 
actions. 

The analytic expressions in Eqs. (l32l) and (1261) for K 
and Ds are profitably used in the following discussion of 
the hydrodynamic function of core-shell particle systems. 


al. found that h d (q) is at all q nearly independent of 
the permeability parameter X x , practically in the com¬ 
plete liquid-phase concentration interval. In extrapolat¬ 
ing their finding to the spherical annulus model as mo¬ 
tivated by our general discussion in Subsec. IHH H(q) 
can be expected to be well represented by 


D. Hydrodynamic function scaling 


Due to the fact that H(q) is given according to Eq. 
(1201) by an equilibrium average over hydrodynamic mo¬ 
bilities, its principal peak location and the wavenumber 
locations of its secondary maxima are nearly coincident 
with those of S(q). As a static equilibrium property, 
S(q) is independent of the hydrodynamic particle struc¬ 
ture and the His in general. This observation has led 
Abade et al. to the following remarkable finding fl 8 j . 
which they analyzed in the context of uniformly perme¬ 
able hard spheres: While the amplitudes of the oscilla¬ 
tions in H(q) are strongly sensitive to the permeability 
(i.e., the hydrodynamic particle structure), the relative 
q-dependence of H{q) is practically permeability inde¬ 
pendent and can be scaled thus to that of no-slip hard 
spheres. To see this quantitatively, consider the so-called 
reduced hydrodymanic function [181 ] . 


hd (q) 


H d (q ) 

\M<1 = 0 )| 5 


(33) 


where H d {q) is the wavenumber-dependent distinct part 
of H(q) introduced in Eq. (1?T1) . 

The reduced hydrodynamic function is defined such 
that h d (q = 0) = —1 and h d (q —I 00 ) = 0. Abade et 


H{q) 


Ds ( 0 , 7 ) 
Dl (7) 


+ K s {q) 


K{<t>, 7 ) 


Ds (</>! 7 ) ' 

Dl ( 7 ) J ’ 


(34) 

where h^ s (q) = h d {q, r ) = 1 ) is the reduced hydrody¬ 
namic function of no-slip hard spheres. The hydrody¬ 
namic particle structure enters into this expression by 
the coefficients K and Ds/Dl only for which we have 
provided accurate analytic expressions. The relative q- 
dependence of H(q) is described by the master function 
h^ s (q) which can be conveniently calculated using the 
semi-analytic <5 7 method for H d (q). The S 7 method has 
been shown, in comparison to Stokesian dynamics simu¬ 
lations, to give quite accurate predictions not only for the 
H d (q) of no-slip hard spheres [43], but likewise for charge- 
stabilized susp ensions with long-range electrostatic in¬ 
teractions [43, [63]. To validate hydrodynamic function 
scaling for the spherical annulus model, in Fig. [FT] we 
present results for h d (q, 7 ) in a broad 7 parameter range, 
obtained using the S'y method in Appendix [B] All curves 
collapse practically on a single master curve which in 
turn nicely agrees with the simulation data for no-slip 
(ns) hard spheres taken from jT§j. 

The sensitivity of H ( q ) on the reduced hydrodynamic 
radius 7 is illustrated in Figure [12] where H(q) has been 
calculated according to Eq. (13H) . The strength of the 
His ceases with decreasing 7 . Notice that H(q) —> 1 for 
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FIG. 12. Hydrodynamic function, H(q), of spherical annu¬ 
lus spheres for <j> = 0.35 and values of y as indicated. The 
solid lines show results by the semi-analytic formula in Eq. 
is. with Ds/Dq and K according to Eqs. (1261) and (1321) . 
respectively, and h^i a (q) calculated using the 5y method. 


7 —> 0 . 



FIG. 13. Reduced hydrodynamic function, hd(q), for a 
Hertzian spheres system of interaction strength e = 10 
(dashed lines) and t = 10 4 (solid lines), respectively, for var¬ 
ious reduced hydrodynamic radii, 7 = a,h/a 3 , as indicated. 
The effective volume fraction is <p B = 0.35. The depicted 
results have been calculated using the S'y method with PY 
structure factor input. Inset: Excerpt of the Hertz potential 
curve for e = 10 (dashed) and e = 10 4 (solid), respectively. 
Wavenumber q and radial distance r are scaled by the effective 
soft diameter, cr s = 2a s , of the Hertz potential. 

To date, the validity of the hydrodynamic function 
scaling was scrutinized for particles with hard-sphere in¬ 
teractions only [13, UK ■ As discussed in Subsec. Ill A1 
the soft Hertz potential in Eq. m is a useful descrip¬ 
tion of the coarse-grained interaction of certain types of 
mechanically soft microgel particles. For this reason, we 


investigate now the scaling of the hydrodynamic func¬ 
tion for suspensions of Hertz particles of different inter¬ 
action strengths e. I 11 Fig. [13] the functions hd(q) of 
Hertz particles are shown for various values of the re¬ 
duced hydrodynamic radius, defined here by 7 = ah/a s 
with a s = a s /2 denoting the effective soft radius of the 
Hertz potential. Two largely distinct interaction (soft¬ 
ness) parameters e = 10 and 10 4 are considered, rep¬ 
resenting highly soft and weakly soft particle systems, 
respectively. The inset depicts the respective shapes of 
the Hertz potential. For e = 10, there is a significant 
likelihood of finding two particles at a distance smaller 
than cr s , as quantified by values of g(r < a s ) significantly 
larger than zero. Even then the HRM remains applica¬ 
ble, provided the hydrodynamic structure of the actual 
soft particles is not significantly distorted away (on aver¬ 
age) from spherical symmetry during particle interpene¬ 
tration. 

The curves for hd{q) in Fig. [13] have been obtained us¬ 
ing the HRM-based Sy method with the structure factor 
input for the Hertzian spheres calculated in PY approx¬ 
imation. The figure shows that hydrodynamic function 
scaling applies to Hertz model particles for a broad soft¬ 
ness range, with the shape of the 7 -independent master 
curve for hd(q) depending on the softness parameter. The 
scaling behavior of H ( q ) can be expected to hold also for 
other soft pair potentials. 


E. High-frequency viscosity 

So far diffusion related properties have been addressed 
only. We discuss next a short-time rheological property, 
namely the high-frequency-limiting suspension viscosity, 
t/qq , which linearly relates the average suspension shear 
stress to the applied rate of strain in a low-amplitude, 
high-frequency shear experiment. Like the short-time 
diffusion properties discussed before, is a quantity of 
purely hydrodymanic origin, influenced by direct particle 
interactions through the equilibrium averaging only. For 
particles acting hydrodynamically like points (7 = 0 ), it 
reduces thus to the shear viscosity, 770 , of the suspending 
Newtonian fluid. 

The high-frequency viscosity should be distinguished 
from the zero-frequency viscosity [80|, 

V (</>) = Voo (</>) + A 77 ((/)) , (35) 

determined in a steady low-shear experiment. The zero- 
frequency viscosity 77 has an additional contribution, 
A 77 > 0, originating from the relaxation of the shear- 
distorted dynamic particle cage formed around each par¬ 
ticle. The shear relaxation part A 77 is influenced both 
by direct and hydrodynamic interactions, with the con¬ 
sequence that for strongly correlated colloidal particles 
it is substantially larger than 7700 . An analytic method 
of calculating the long-time transport property 77 is pre¬ 
sented in Subsec. IV Bl 
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FIG. 14. High-frequency viscosity, 7 of spherical annulus 
particles in dependence on 0, for values of 7 as indicated. 
Filled symbols: Simulation data taken from [ 2 ll . :67|. Solid 
lines: Generalized Saito formula in Eqs. (CISll - (IMD- Dashed 
lines: PA scheme results with VW-PY input for g(r). Dotted 
lines: 5 7 scheme results with VW-PY input for S(q). 


In [H, 0 , a generalized Saito formula for the high- 
frequency viscosity of permeable spheres with hard-core 
interactions has been introduced which for 0 < 0.5 gives 
results in good agreement with simulation data. In the 
framework of the spherical annulus model, the formula 
reads 


Voo ( 4 >, 7 ) 
Vo 


1 + [77] ( 7 )<t>~ 
1 


l + g(0,7) _ 

I [v] (1 + s ( 0 , 7 )) 


(36) 

It expresses Voo in terms of the intrinsic viscosity, 
[v\ ( 7 ) = (5/2) 7 , depending on the hydrodynamic par¬ 
ticle structure only, and the Saito function S (0, 7 ). The 
latter is approximated linearly in 0 as 


s (^ 7 ) 


( V(7) 
VW(7) 


: W(7) 


(37) 


where Ay ( 7 ) is the first-order virial coefficient in the ex¬ 
pansion of Voo/Vo in powers of 0. Numerical values for 
the first virial coefficient of spherical annulus particles 
are given in (22). For 7 > 2/3, these values are well 
represented by the polynomial 

A v ( 7 ) = 5.0021 - 39.2797+ 143.1797 2 (38) 

- 288.2027 3 + 254.5817 4 . 

In Fig. 1141 the predictions for r)oo by the generalized 
Saito formula in Eqs. (1551) - (1551) are compared with 
existing simulation results [ 2 l| for the spherical annulus 
model. There is good agreement with the simulation in 
the displayed liquid-phase concentration range. For fixed 
concentration n and fixed hard-core radius a, the viscos¬ 
ity increases with increasing ah, he. increasing 7 , owing 


to the enlarged dissipation. Additionally shown in the 
figure are results for 77 ^ by the PA and 67 methods de¬ 
scribed in Appendices lAl and iBl respectively. Like in the 
PA scheme for short-time diffusion properties, two-body 
His contributions to Voo are fully accounted for but three- 
body and higher order contributions have been neglected. 
The PA scheme is in good agreement with the simulation 
data for 0 < 0.2 only. For small 7 , the applicability of 
the PA method extends to somewhat larger 0. We at¬ 
tribute this first to the weaker hydrodynamic interactions 
for 7 < 1 , and second to the fast 0 (l/r 6 ) asymptotic de¬ 
cay of the shear mobility function associated with 7700 
(see Appendix lAl). While the 67 scheme viscosity predic¬ 
tions for no-slip spheres are in better agreement with the 
simulation data than those by the PA scheme, for 7 < 1 
the ($7 scheme consistently overestimates the high- fre¬ 
quency viscosity. Quite interestingly, for all 0 < 0.5 the 
dy prediction for 7700 scales to good accuracy with the 
hydrodynamic volume fraction 0 h = 7 3 0 . 



FIG. 15. Test of the short-time GSE relation in Eq. (1391) 
for spherical annulus particles with values of 7 as indicated. 
Solid lines: Product function Asv of the generalized Saito 
expression in Eq. (l36l) for 700/70 and Eq. (l26l) for Ds/Dq. 

As straightforward applications of the generalized 
Saito formula and Eq. (l26l) describing 7 ^ (</>, 7 ) and 
Ds{4>, A) in the spherical annulus model, we analyze next 
the validity of the short-time generalized Stokes-Einstein 
(GSE) relation, 


AsvW>,7)~1, (39) 


with the short-time GSE function 


Asv(</ ; 7 ) 


Ds{<t>,l) 7700 ( 0 , 7 ) 

^ 0 ( 7 ) Vo 


(40) 


Eq. (l39l) expresses that Ds(4>, 7 ) should be propor¬ 
tional for all concentrations to the inverse of Voo (4>, 7 ) ■ 
This relation is trivially fulfilled at infinite dilution 
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where it reduces to the single-sphere translational Stokes- 
Einstein relation for a hydrodynamically structured col¬ 
loidal sphere. The approximate validity of this relation 
would be quite useful, since rjoo can then be determined 
more easily, and using a smaller amount of particles, 
by a dynamic scattering experiment instead of a rheo- 
mechanical measurement. This is why GSE relations in¬ 
cluding the present one have been thoroughly subjected 
to theoretical explorations, for particulate systems in¬ 
cluding permeable hard spheres [ 20 ] and charge-stabilized 
particles [H, [67| . 

An exact GSE relation is reflected in Fig. (fT5l) by a hor¬ 
izontal straight line of unit height. However, for all con¬ 
sidered values of 7 , significant deviations from Agy = 1 
are observed at larger volume fractions. The deviations 
are largest for no-slip hard-core particles where the His 
are strongest, amounting to about 40% at cj> = 0.45. For 
concentrations (j) < 0.4, the displayed curves for As v(4>) 
are nearly straight lines, characterized by the linear co¬ 
efficient, Asv( 7 ), in the expansion of A sv to linear order 
in </>. The linear concentration coefficient derived from 
our analytic expressions for D$ and rj oo is given by the 
polynomial 

Asv( 7 ) = 0.6685 + O. 32 OI 7 + O (y 2 ) . (41) 

For a given hydrodynamic particle model, the values 
of 7 which should be used in calculating Ds and r?oo, 
respectively, are actually different if the O(L^j) 2 correc¬ 
tions to a,hj in Eq. (ThU) cannot be neglected. However, 
this does not affect our general conclusion that the GSE 
relation is violated, as illustrated by the curves in Fig. 
(1151) where for simplicity the same 7 values were used in 
calculating Ds and ijoo- 

According to Fig. (USD, the ordering relation Asv(^) > 
1 is obeyed by particles with pure hard-core interac¬ 
tions. As shown in [67| . the same ordering applies to 
charge-stabilized suspensions. In contrast, the long-time 
product function Alv(</>) = (Dl/Dq) x {q/r/o) relating 
long-time self-diffusion coefficient Dl to zero-frequency 
viscosity q has been shown for no-slip hard-sphere and 
charge-stabilized suspensions to fulfill the opposite or¬ 
dering Alv(^) < 1 for </> > 0 [H]. 


IV. THEORY VERSUS EXPERIMENT 

We demonstrate now the accuracy of our easy-to-apply 
toolbox methods of calculating short-time diffusion prop¬ 
erties by analyzing DLS measurements by Eckert and 
Richtering [43] on non-ionic PNiPAM microgels in DMF. 
As discussed in Subsec. Ill Af the microgel particles be¬ 
have as hard spheres as far as their static properties are 
concerned. On modeling the microgels also hydrodynam¬ 
ically as no-slip hard spheres with ah = a , and on basis of 
bare 5 7 method results for H(q), Eckert and Richtering 
came to the conclusion that short-time particle diffusion 
is underestimated by the no-slip hard-sphere model. This 


suggests that the non-uniform crosslinker-density of the 
microgels should have a significant hydrodynamic effect. 

To account for this effect, we model here the microgels 
as spherical annulus particles, and determine H{q) us¬ 
ing the scaling Eq. (1341) in conjunction with the analytic 
expressions in Eqs. (El and (® for K and Ds, re¬ 
spectively. In the ($7 method calculation of h(q) which 
enters into the scaling expression of H(q), we use the 
HS-PY structure factors depicted in Fig. [l] with the 
hard-core radius a = 120 nm. The only adjustable pa¬ 
rameter in our model is thus the reduced hydrodynamic 
radius 7 = ah / a . 



FIG. 16. Experimentally deduced hydrodynamic function 
of PNiPAM microgels in DMF (filled symbols, taken from 
0) compared with the theoretical predictions (solid lines) 
for the spherical annulus model, using 7 = 0.97 and Eq. El 
for H(q) combined with Eq. (1321) for K and Eq. (12611 for 
Ds- Dashed line: Theoretical prediction for non-permeable 
particles (7 = 1) at <j> = 0.5. The wavenumbers are scaled by 
the hard-core diamter a = 240 nm. 

In Fig. m ou r theoretical results for H(q) are pre¬ 
sented and compared with the experimental findings in 
[0 . The latter have been obtained indirectly from divid¬ 
ing the DLS first cumulant data for the diffusion function 
D(q) by the VW-PY S(q) taken at q = 0 (c.f. Eq. (flHll b 
Using the constant ratio 7 = 0.97, we obtain very good 
agreement between theory and experiment for all volume 
fractions. Our finding of a concentration-independent hy¬ 
drodynamic radius ah = 0.97 x a points to the consis¬ 
tency of our analytic method of calculating H(q), since 
as an intrinsic particle property, ah should not depend 
significantly on the volume fraction. While this holds 
for the strongly cross-linked non-ionic microgels consid¬ 
ered here, for weakly cross-linked ionic microgels in the 
swollen-state temperature range a significant size shrink¬ 
age with increasing concentration is observed [lj. 

The deduced hydrodynamic microgel radius is only 3% 
smaller than the excluded volume radius, correspond¬ 
ing to a likewise small value, \ x = 0.029, of the re¬ 
duced fluid penetration length. This exemplifies the com- 












16 


mon experimental situation of ah ~ dhj, with a rela¬ 
tive correction to the flat plane value ahj being here of 

O ~ 10 " 3 small. 

The small microgel permeability nonetheless signifi¬ 
cantly affects H(q), in particular at larger concentrations. 
This is shown in Fig. [16] for (j> = 0.5 by the comparison 
with the hydrodynamic function for zero permeability 
(dashed curve): The residual particle permeability en¬ 
larges the sedimentation velocity by more than 100 %, 
and the self-diffusion coefficient by more than 30%. 



FIG. 17. DLS data for the reduced diffusion function D(q)/D q 
(filled symbols) taken from 23] j in comparison with our theo¬ 
retical predictions for 7 = 0.97, obtained from the multiplying 
the sperical annulus results for H(q) depicted in Fig. 1161 bv 
the VW-PY values of S{q = 0). 


In Fig. [TT] our theoretical results for H(g)/Dg are 
plotted together with the DLS data 23 for the same 
quantity. The theoretical curves have been obtained from 
dividing the spherical annulus H{q )'s depicted in Fig. [16] 
by the associated VW-PY structure factors, S(q), of hard 
spheres in accordance with Eq. (Il9l) . At q = 0 , the VW- 
PY S(q) reduces to the Carnahan-Starling expression for 
the reduced osmotic compressibility, 


S(q = 0) 


(1 + 2(f)) 2 + <f> 3 (cf> ~ 4) ’ 


(42) 


valid in the full fluid-phase volume fraction range of hard 
spheres. The agreement between theoretical and experi¬ 
mental diffusion functions is very good in the intermedi¬ 
ate wavenumber range including the principal peak po¬ 
sition, q mi of S(q) where D(q) is minimal, and also for 
larger wavenumbers. At small q values and large volume 
fractions, the experimental data are overestimated. Even 
considering the remaining small -5 deviations, the here 
presented theoretical results for D(q) are in distinctly 
better agreement with the experimental data than the 
earlier ones presented in 23 where permeability effects 
were not included. 


The deviations in D(q) at small -5 can be partially at¬ 
tributed to the high sensitivity of the inverse compress¬ 
ibility factor, 1/5(0), multiplying H( 0) in Eq. (1^1) . on 
the residual softness of the microgels. This is demon- 



FIG. 18. Influence of particle softness (elasticity) on the in¬ 
verse compressibility factor 1/5(0), for hard-sphere and Hertz 
potential particles, plotted as a function of <j>. The Carnahan- 
Starling (CS) (black solid line) and PY (red solid line) results 
for hard spheres are compared with the PY-based pedictions 
for Hertz model particles with softness parameters e = 10 4 
(dashed line), 10 2 (dashed-dotted line), and 10 (dotted line). 


strated in Fig. [TS] where the concentration dependence of 
1/5(0) for hard spheres is compared to that of the Hertz 
potential system for the strongly distinct softness param¬ 
eters e = 10 4 , 10 2 and 10. According to the figure, even 
a small residual softness characterized by e = 10 4 sig¬ 
nificantly enlarges the osmotic compressibility for large 
volume fractions, with 1/5(0) being lowered accordingly. 
Viewed on the scale of the structure factors in Fig. [2] 
the small -5 differences are not resolved since 5(0) is very 
small for large concentrations. While the smaller factor 
1/5(0) in the Hertz model would improve the agreement 
with the experimental D{q) at small q , we recall in re¬ 
ferring to Fig. [2] that a somewhat larger volume fraction 
cf> s than for hard spheres is required in the Hertz model 
for an equally good fit of the experimental S(q). If the 
enlarged volume fraction is accounted for, there remains 
a small final reduction in 1/5(0) only. An additional 
cause for the small -5 differences can be size polydisper- 
sity. It was shown in [ 8 l|, that a small degree of polidis- 
persity significantly enlarges in concentrated suspensions 
the measured diffusion function at small 5 . 

In summarizing, we conclude that our HRM based 
toolbox methods reproduce the short-time diffusion prop¬ 
erties of the non-ionic microgel suspensions very well, and 
with little numerical effort. The non-homogeneous cross¬ 
linker density is accounted for by a hydrodynamic radius 
which is only three percent smaller than the excluded 
volume radius. 
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V. LONG-TIME TRANSPORT PROPERTIES 


A. Long-time self-diffusion 


Long-time colloidal transport properties such as Dl 
and 77 characterize the particle dynamics on time scales 
t tjj during which the particle configuration has 
changed significantly. Different from short-time proper¬ 
ties, they are influenced by thermally driven microstruc- 
tural relaxations depending on direct and hydrodynamic 
interactions alike. This renders a first-principles calcula¬ 
tion of long-time properties demanding, both in theory 
and simulations, in particular when the salient His are 
accounted for. Consequently, in most simulation studies 
of the concentration dependence of Dl and n, the influ¬ 
ence of the His has been ignored [82| [HI j84[ • The few 
existing three-dimensional simulation studies where His 
are included have been focused on Brownian hard spheres 
with no-slip BC [85M87j ] . Therefore, a theoretical scheme 
is on demand allowing for the approximate calculation of 
Dl and 77 for suspensions of hydrodynamically structured 
particles. 

In the following, we present such a scheme which as 
a bonus requires only little numerical effort. It is based 
on the HRM and a factorization approximation method 
proposed originally by Medina-Noyola for self-diffusion 
[351 1 • We point out that the HRM error estimation for 
short-time properties in Eq. (l25l) applies also to long¬ 
time properties including Dl/Dq and 77 , and also the 
dynamic structure factor S(q,t) for arbitrary correlation 
times. This follows from general expressions for long¬ 
time transport coefficients and S(q,t) which have been 
obtained using the Mori-Zwanzig projection operator for¬ 
malism in conjunction with the many-particle generalized 
Smoluchowski equation for the configurational probabil¬ 
ity distribution function a The crucial fact to notice 
here is that the hydrodynamic mobility tensors entering 
into the GSmE have no explicit time dependence. 

We exemplify the error estimate for Dl by starting 
from the configurational average expression [3l|, 

D l = D s + AD, (43) 

with 

AD = -<(q-V 1 )Os 1 (q-V 1 )) . (44) 

The slowing effect on Dl by the dynamically restructur¬ 
ing cage of next neighbors formed around each particle 
is embodied in the negative valued relaxation contribu¬ 
tion, AD, to Dl implying Dl < Ds for cj> > 0. In 
Eq. (Ell . Vi = ObRi is the coarse-grained velocity of a 
representative particle 1 whose Brownian motion is con¬ 
sidered, q = q /q, and Og(X) the so-called backward 
Smoluchowksi differential operator generating the time 
evolution of the probability distribution function. The 
operator inverse is denoted by Og 1 . The only informa¬ 
tion about Ob needed here is its linear dependence on 
the mobility tensors Hij. From this and Eq. (Ell , it fol¬ 
lows that the error introduced on approximating Dl / Dq 

by Dl-krm/DI is of O ((££,/) 2 )- 


Using Eq. m, Dl can be written as 


Dl(4) 


Ds (</>) 


A dm; 

DsW _ 


(45) 


where in the term in brackets, the explicit dependence on 
Dg has been scaled out. According to arguments first put 
forward by Medina-Noyola [35|, and subsequently sub¬ 
stantially elaborated by Brady also regarding the zero- 
frequency viscosity [36l l37j] , the factor function in brack¬ 
ets is not only scale invariant with respect to Ds, but 
for hard spheres it is to decent approximation also inde¬ 
pendent of the His. This implies the no-HI factorization 
approximation, 


1 + 


AD 

~Dj 


Dl 

Dl 


-HI 


(46) 


where the bracket term is approximated by a purely 
structural factor determined by excluded volume inter¬ 
actions only. For known Dg, the problem of calculating 
Dl is thus simplified to the calculation of the reduced 
long-time self-diffusion coefficient without His. The de¬ 
pendence of Dl on the His, and the hydrodynamic par¬ 
ticle structure, is embodied here in Dg alone. 

Brownian dynamics simulation results for 
[Di(^)/Dg]no-Hi by Hinsen and Cichocki [Hj and 
Moriguchi [ 88 | are depicted in Fig. [19] In the fluid-phase 
concentration regime <f> < <fif, where cj>f = 0.494 is the 
volume fraction at freezing, the simulation data are well 
described by the polynomial least-square fit, 


= 1 20 -(- 1 . 2720 2 — 1 . 9510 3 ( 47 ) 

no-HI 

where the exact first-order virial coefficient, Xl = 2 , for 
a hard-sphere suspension without His has been incorpo¬ 
rated. The figure depicts furthermore the analytic ap¬ 
proximation, 



D l 


n* 
u 0 


o-HI 


1 

1 + 2 <t>g(a+■,<!>) ’ 


(48) 


given by Brady (36}, where the structural factor is ex¬ 
pressed in terms of the RDF contact value, g(a + ',4>), of 
hard spheres. The contact value for the fluid-phase con¬ 
centration range is to high accuracy described by the 
Carnahan-Starling expression [89| 


g(a + ;4>) = - -(49) 

( 1 - 0) 3 

with g(a + \(j)f) = 5.81. Eq. (H51) incorporates the exact 
first-order virial coefficient, A l = —2. Moreover, near 
random closed packing at <p lcp « 0.64 where a metastable 
hard-sphere fluid gets jammed, on basis of results for 
g(a + ;<f) by Rintoul and Torquato i89j it predicts that 
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[-DzJno-Hi diminishes linearly like 0.59 x (1 — (j)/(j ) TCp ). 
Since Dg vanishes likewise linearly in case of no-slip 
hard spheres, the quadratic scaling prediction Dl ~ 
(1 — <f>/(j> rC p ) 2 near random closed packing is obtained 
in the factorization approximation. 



FIG. 19. Reduced long-time self-diffusion coefficient, 
\Dl/DQ , of colloidal hard spheres without His. Filled 
circles and squares: Brownian dynamics simulation results 
by Hinsen and Cichocki f§2} and Moriguchi [88|, respectively. 
Solid line: Polynomial fit in Eq. m- Dashed line: Eq. (1481) 
with Carnahan-Starling contact value input. 

The performance of the no-HI factorization approx¬ 
imation for the Dl of no-slip spheres on basis of Eq. 
J26| for Ds at 7 = 1 and Eq. ([47} for [D l /D^\ no _ H i> 
is documented in Fig. [20] by the comparison with DLS 
data by van Megen et al. J90l H3 , and simulation re¬ 
sults by Phung et al. [85| with His included. The 
Stokesian dynamics simulation data for Brownian hard 
spheres by Phung et al. have been obtained for a small 
number (N = 27) of particles, and for a small albeit 
non-zero shear Peclet number Pe = 0.01. The overall 
agreement with the experimental and simulation data 
is quite good. The factorization approximation gives 
A l = —3.831 for the first-order virial coefficient of hy- 
clrodynamically interacting no-slip spheres, while its cor¬ 
rect numerical value is given by A l = —2.1. The initial 
low-concentration decrease of Dl is thus overestimated. 

The good performance of the no-HI factorization ap¬ 
proximation for no-slip hard spheres gives support to 
its straightforward extension to hydrodynamically struc¬ 
tured particles, by using for Dg/Dg now the analytic 
expression in Eq. (1261) for spherical annulus spheres. 
Our results for Dg/D\ based on this extended factor¬ 
ization scheme are shown in Fig. [5TJ With decreas¬ 
ing 7 , the slowing down effect of the His on Dl di¬ 
minishes. In the limit 7 —>• 0 of hard spheres act¬ 
ing hydrodynamically as point particles, the long-time 
self-diffusion coefficient in the absence of His is recov¬ 
ered. According to Cichocki and Fclderhof [7dj |. the lin¬ 
ear concentration coefficient of Dl/Dq for small 7 is 



FIG. 20. Reduced long-time self-diffusion coefficient, Dl/Dq, 
of no-slip colloidal hard spheres. Filled circles: Experimental 
data by van Megen et a l. I90L lol ]. Filled squares: Simulation 
data by Phung et al. [ 83 . Solid line: No-HI factorization 
approximation in Eqs. with the short-time factor 

Ds/Dq according to Eq. (1261) for 7 = 1. 



FIG. 21. Generic influence of the intra-particle hydrodynamic 
structure on Dl/Dq, quantified in the no-HI factorization ap¬ 
proximation, Eqs. G3) and 03, using for Ds/D q the spher¬ 
ical annulus model formula in Eqs. (HO) and (|2T|). Several 
values of 7 are considered as indicated. 


A 1 ,( 7 ) = —2 [l — 1.0317 + O.III 7 2 + C ) ( 7 3 )], where the 
linear and quadratic terms in 7 are due to the so-called 
Oseen long-distance His contribution to the relaxation 
part AD. The contributions of the short-time part, Dg, 
to Dg appear first in quadratic order in 7 . While this de¬ 
scribes quantitatively how Dl approaches [-Dl] n o-Hi for 
small 7 , we reemphasize that 7 > 0.8 for most hydrody¬ 
namically structured colloidal particles. 

The no-HI factorization approximation predicts the ra¬ 
tio, .Dl/Ds, of long-time and short-time coefficients to 
be independent of His and hydrodynamic particle struc- 


























ture, with value equal to [Dl/ .Do]no-Hi- This implies with 
Eqs. fl5l) and (PTO!) that 
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(^)(0/)« 01 , (50) 

in good accord with the Lowen-Palberg freezing criterion 
value of about 0.1. Thus, a universal freezing value of 0.1 
is predicted not only for colloidal suspensions with differ¬ 
ent pair potentials, but also with different hydrodynamic 
intra-particle structures. For pair interactions character¬ 
ized by a single length scale, the dynamic Lowen-Palberg 
criterion has been shown to be equivalent to the static 
Hansen-Verlet freezing criterion for the value S(q m ) of 
the structure factor peak height. 



B. Zero-frequency viscosity 


Analogous to Eq. (TT5l) for Dl, in Eq. (l35l) for the low- 
shear zero-frequency viscosity 7?, we factor out the high- 
frequency (short-time) contribution r/oo according to 


v{ 4 >) 


Voo ( 0 ) 


' At? (0) ~ 

. Voo(<P)\ ’ 


(51) 


with the term in brackets expected to be approximately 
independent of the His. This suggests the no-HI factor¬ 
ization approximation, 


A?? (0) 

Voo (0) 


At? (0) 


L Vo 


no-HI 


(52) 


where [AT?] no _Hi is the shear relaxation viscosity part 
without His. In this approximation, the His are assumed 
to affect 77 only by means of the factored out t?oo in Eq. 
(EH). The neglect of His simplifies the calculation of the 
shear relaxation viscosity part considerably. Following 
works by Brady [3f| [37|, an analytic estimate of At? for 
no-slip hard spheres without His is given by 


At? 

Vo 


^(/) 2 g(a + ;cl)) 


(53) 


- no-HI 


This estimate combines the exact low concentration 
limit, 2.4 0 2 + 0(0 3 ), of [A7?] no -Hi with its divergence 
at random closed packing according to [AT?] no _Hi ~ 
(1 — 0 / 0 rc p) -1 , triggered by the divergence of the hard- 
sphere contact value. Together with the likewise linear 
divergence of t?oo for no-slip hard spheres, a quadratic di¬ 
vergence 77 ~ (1 — 0 / 0 rcp ) -2 is thus predicted for the 
zero-frequency viscosity. For hydrodynamically struc¬ 
tured particles where ah < a , there are no diverging lubri¬ 
cation forces for spheres in contact. The high-frequency 
viscosity remains then finite at random closed packing, 
and the particles can still rotate individually. Moreover, 
t? diverges only linearly as rj(4>, 7 < 1 ) ~ (1 — 0 / 0 rcp ) . 

In Fig. [221 the outcome of Eq. (1551) for [AT?] no _ H i is 
compared to Brownian dynamics simulation results with¬ 
out His by Foss and Brady [H], [H]. Up to 0 « 0.35, the 


FIG. 22. Reduced shear relaxation viscosity part without 
His, [At? / t?o] no-Hi , of a suspension of Brownian hard spheres 
in dependence of 0. Filled circles: Brownian dynamics (BD) 
simulation data by Foss and Brady (84j[. Solid line: Analytic 
estimate in Eq. (1531) . Dashed line: Semi-empirical fit in Eq. 

EH. 


simulation data are decently well represented by the an¬ 
alytic expression, but the steep rise of [AT?] no _Hi at large 
volume fractions is not reproduced. The simulation data 
are well captured for all 0 by the semi-empirical expres¬ 
sion 


A?? 

W (1 - 

7.085^ + 20.1S20 2 ) 

. Vo . 

no-HI 

t 1 -^) 



(54) 


combining the exact quadratic-order concentration de¬ 
pendence with the linear order divergence at random 
closed packing. Throughout this work, we restrict 
our analysis to the equilibrium fluid-phase concentra¬ 
tion regime 0 < 0.5, while the viscosity simulations in 
[83l HU, HH have been extended to the metastable fluid 
concentration regime 0f < 0 < 0 rcp where crystallization 
is kinetically suppressed. Regarding the shear relaxation 
viscosity part with His included, Brady has proposed the 
following approximate scaling expression [36], [37[, 


At? _ 12 2 g(o~+; 0) 
7700 5 ^ A sv (0) 


(55) 


where the influence of the His on the ratio Arj/rjoo is 
solely embodied in the short-time GSE function Asv de¬ 
fined in Eq. m . For no-slip hard spheres, Asv is well 
represented for 0 < 0.4 by Asv ~ 1 + 0.67 <j> according to 
Eq. (PHI) , while without His Asv is equal to one. Con¬ 
sequently, Eq. (l55l) predicts Arj/r/oo to be only mildly 
affected by the His, giving some credit to the no-HI fac¬ 
torization approximation in Eq. (1521) . Eq. (1551) was 
obtained from arguing that the adequate diffusion time 
scale in a concentrated suspension is a\/D s instead of 
a h/D 0 , and from using a low-concentration estimate of 
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FIG. 23. Zero-frequency viscosity, 77 / 770 , of no-slip Brownian 
hard spheres with His. Solid and dashed black lines: No- 
HI factorization predictions using ?7oo(<T>,7 = 1) according to 
Eq. (1361) . and [?7/?7o]no-Hi according to Eqs. (1531) and EU) , 
respectively. Solid and dashed red lines: Brady’s (modified) 
scaling approximations using Eqs. (1551) and (1561) for Ag/rioo, 
respectively, and Eq. (13611 for the factored out r/oo. The factor 
1/Asv in the (modified) scaling approximation is calculated 
using Eq. (HO) for r/oo and Eq. (1261) for Ds■ Dotted line: 
Eq. m- Filled symbols: Experimental data by Segre et al. 
and Weiss et al. [92. [9311 . Open symbols: Stokesian dynamics 
(SD) simulation results for Brownian hard spheres by Foss 
and Brady [86[ and Phung [85} . 


the weakly shear-distorted stationary pair distribution 
function with the prefactor g(a + ;(f>) preserved [36], A 
detailed discussion of the approximatio ns g oing into Eq. 
(1521) is given by Lionberger and Russel [94 195|. 

On using in place of Eq. (l53l) the semi-empirical fitting 
expression in Eq. (1541) as the non-hydrodynamic factor 
in Eq. (l55l) , Brady’s scaling relation is modified to 


Ag _ f </- 2 (1 - 7.085</ + 2O.1820 2 ) 


(56) 


The two here considered variants of the no-HI scal¬ 
ing approximation of 77 consist of using Eqs. ED and 
m, respectively, as input for the ratio Ag/goo in the 
bracket term in Eq. ED, in conjunction with the accu¬ 
rate generalized Saito formula in Eq. (1331) used for r^. 
In addition, Brady’s scaling expression for 77 and its mod¬ 
ification consist of approximating Ag/rjoo in the bracket 
term in Eq. ED by Eqs. (1551) and (1561) . respectively, with 
the generalized Saito formula used again for the factored 
out high-frequency viscosity. The hydrodynamic factor 
1/Asv in the two scaling expressions is calculated ana¬ 
lytically using Eq. (1531) for Ds , and the generalized Saito 
formula for 7]^. 

The results for 77 (^ 1,7 = 1) by the four inter-related 
analytic approximations are depicted in Fig. 1231 They 
are compared with experimental data by Segre et al. |92[ 


and Weiss et al. [93(, and Stokesian dynamics (SD) simu¬ 
lation data for Brownian hard spheres with His included 
by Foss and Brady [86} and Phung at al. [§5}. The lat¬ 
ter have been obtained for a small number (TV = 27) of 
particles in the basic simulation box. In [8f|, A77 was de¬ 
duced using a general Green-Kubo formula for the shear 
stress correlation function of hydrodynamically interact¬ 
ing particles [961 ] . The viscosity curves by all four approx¬ 
imations compare overall quite well with the simulation 
data. Depending on the used approximation for A 77 / 7700 , 
deviations are visible at intermediate and large volume 
fractions. The neglect of His in Arj/rioo results in an over¬ 
estimation of the simulations data at large <f>, but the 
large-</> experimental viscosity data are well described. 
On considering the simulation data to be more trustwor¬ 
thy than the experimental data, owing to experimental 
polydispersity effects and difficulties in determining the 
precise volume fraction, the modified scaling expression 
by Brady in Eq. (1561) provides the overall best description 
of the SD data for 77, slightly better than Brady’s original 
scaling expression in Eq. (l55l) . We use the modified scal¬ 
ing expression in the following discussion of hydrodynam¬ 
ically structured particles. Note that the second-order in 
concentration coefficient of 77 is predicted by all four fac¬ 
torization approximation variants as 5.01 + 2.4 = 7.41, 
while the exact coefficient is equal to 5.931 (97}. This 
low -(f> difference is not resolved on the scale of Fig. [53] 
The figure includes additionally the viscosity predic¬ 
tion by the formula 


!- 0-40 + 0 . 222 / 

vW/Vo 2 ' (°0 

(1-0) 

with (j> = (f>/(f> IC p , which incorporates the first two known 
virial coefficients in 77/770 = 1 -|- 2.5</>-F 5.91</> 2 -+-Cf7(</) 3 ) and 
a quadratic divergence of 77 at </> rcp . As seen in Fig. [ 53 ] it 
describes the simulation and experimental data well for 
(j) < 0.35, but it strongly underestimates them at larger 
</>■ 

Akin to long-time self-diffusion, we can straightfor¬ 
wardly extend our analysis to hydrodynamically struc¬ 
tured particles, by using the generalized Saito expression 
in Eq. (1361) for r]oo((f>, 7), in combination with the mod¬ 
ified Brady scaling expression in Eq. (1561) for Arj/rjoo. 
The viscosity predictions for different reduced hydody- 
namic radii 7 and volume fractions extending up to <f>{ 
are depicted in Fig. [24] Note the pronounced reduction 
of 7) with decreasing 7, owing to the reduced dissi¬ 
pation, inherited from the similar behavior of 77oo(</>, 7). 
In the limit 7 —>- 0, the zero-frequency viscosity reduces 
to rj((f >, 0) = 77o + [A77] no _Hi. In spite of being approximate, 
the analytic modified scaling expression for rj((f), 7) in Eq. 
(1531) can be expected to be useful for a quick analysis of 
experimental viscosity data of hydrodynamically struc¬ 
tured colloidal suspensions. For the PNiPAM microgels 
in DMF, e.g., where 7 = 0.97 has been deduced, a signif¬ 
icant reduction both in 77 and rjoo is predicted relative to 
the corresponding viscosities of non-permeable particles. 
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FIG. 24. Predictions for the reduced zero-frequency vis¬ 
cosity, rf/r/o, of hydrodynamically structured particles, 
based on the modified Brady scaling expression for A 77/7700 in 
Eq. (1561) . and the generalized Saito formula in Eq. (1361) for 
7 ). The hydrodynamic factor Asv(</>, 7) is calculated 
using Eq. (1551) for UnofaW). and Eq. (1551) for flsfd.'yl/flnhl. 
Several values of 7 are considered as indicated. 

It will be interesting to compare our zero-frequency and 
high-frequency viscosity predictions with future viscosity 
measurements on non-ionic PNiPAM in DMF systems. 


VI. CONCLUDING REMARKS 

We have presented a toolbox of methods for calcu¬ 
lating short-time and long-time transport properties of 
suspensions of spherical particles with intrinsic hydrody¬ 
namic structure. The analytic scaling expressions given 
in the paper combine high accuracy predictions with a 
facile implementation, and they apply to the full liquid- 
phase concentration regime. They are useful to exper¬ 
imenters for a fast yet precise data analysis of scatter¬ 
ing and rheo-mechanical experiments. We have high¬ 
lighted this in our analysis of SLS and DLS experi¬ 
ments on non-ionic PNiPAM microgels in DMF. By a 
detailed comparison with Hertz potential calculations, 
we have shown that the microgels behave statically to 
good accuracy as hard spheres, and hydrodynamically 
as permeable spheres with a reduced penetration length 
X x ft 0.03 corresponding to a reduced hydrodynamic ra¬ 
dius of 7 = 0.97. It will be interesting to scrutinize 
the analytic viscosity expressions for r ] 00 and rj of hy¬ 
drodynamically structured, rigid spheres against future 
rheo-mechanic measurements on non-ionic PNiPAM mi¬ 
crogels. 

The analytic toolbox expressions for the collective dif¬ 
fusion coefficient and the zero-frequency viscosity have 
been already profitably used in a recent crossflow ultra- 
filration study of solvent-permeable nanoparticles sus¬ 
pensions EH- The short-time transport coefficient ex¬ 


pressions can serve also as input in the calculation of 
frequency- and time-dependent transport properties on 
basis of mode-coupling theory and dynamic density func¬ 
tional theory methods were His are included [98, 99j. 
They are likewise useful as input in the context of the em¬ 
pirically observed time-wavenumber factorization scaling 
of Sjq , t) at q ~ q m for intermediate to long correlation 
times [ 10 d 1 10 lT | . 

The intra-particle hydrodynamic structure has been 
accounted for using the hydrodynamic radius model 
where the particles are described hydrodynamically as 
no-slip spheres, characterized by a hydrodynamic radius 
derived from a single-particle transport property for un¬ 
changed direct interactions. In spite of its simplicity, the 
HRM is universally applicable since corrections terms are 
usually quite small, i.e. of quadratic order in the reduced 
slip length. In many particulate systems, these correc¬ 
tions are negligible, and a unique hydrodynamic radius 
dhj can be used independent of the considered transport 
property. As we have shown in comparison with exist¬ 
ing computer simulation results, the HRM based scaling 
expressions for hard-sphere-like particles are decent ap¬ 
proximations also for strongly structured particles char¬ 
acterized by values of 7 significantly smaller than one, 
provided ah is deduced from an associated single-particle 
transport coefficient, namely from [77] in viscosity calcula¬ 
tions, and Dq in self-diffusion and sedimentation calcula¬ 
tions for concentrated suspensions. While our focus has 
been on monodisperse particle systems, it is feasible to 
generalize the HRM to size polydisperse particles using 
an appropriate particle size distribution. In future work, 
it will be rewarding to search for possible extensions of 
the scaling relations to polydisperse colloidal systems and 
mixtures. 

Most results presented in this work are for hydrody¬ 
namically structured, stiff particles with pure hard-core 
interactions where the HRM reduces to the spherical an¬ 
nulus model. Part of our toolbox methods are likewise 
applicable, with appropriate modifications, to spherical 
particles with short-ranged soft interactions. We have il¬ 
lustrated this in Fig. [El] for the reduced hydrodynamic 
function, hd(q), of soft Hertz potential particles with 
e = 10. According to our calculations based on the Sj 
method, hydrodynamic function scaling remains valid for 
soft spherical particles. Note here that the self-part cor¬ 
rected S"f scheme, and the PA approximation for HRM 
particles, can be used for arbitrary pair potentials. 

We finally point to the invalidity of the presented scal¬ 
ing relations for D$, K, and Dl for particles with long- 
range soft repulsion that is not adequately described by 
an effective excluded volume diameter. Examples in case 
are low salinity suspensions of charge-stabilized particles 
[75} . and quasi-two dimensional systems of m agne tically 
repelling particles at a liquid-gas interface |l 02 | . As 
shown both in computer simulations and experiments, 
Ds and K in these systems follow a fractal ~ (j> 4 ^ 3 and 
~ (f) 1 / 3 concentration dependence, respectively, of zero 
and negative infinite initial slope [63, |75j . Moreover, the 
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factorization scaling relation in Eqs. (l45l) and (l46l) pre¬ 
dicting that D l < [Di] no -Hi does not apply to these 
long-range repulsive systems where, on the contrary, a 
hydrodynamic enhancement of long-tim e sel f-diffusion is 
observed in theory [98[, simulation (Yl ll03| and exper¬ 
iments lM- The short-time self-diffusion coefficient 
in these systems is only weakly affected hydrodynami- 
cally at lower concentrations so that Dg ~ Dq. The 
major effect originates instead from the relaxation part 
AD whose magnitude is lowered by the here dominat¬ 
ing far-field part of the His [98|. While the scaling ex¬ 
pressions are invalid, the hydrodynamic structures of the 
low-salinity charge- stabilized and magnetic particles is 
still well described by the HRM. On basis of the HRM 
the short-time properties of hydrodynamically structured 
charge-stabilized particles can be calculated using the 
self-part corrected <5y scheme in conjunction with the PA 
method for the self-part. 
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Appendix A: Pairwise-additivity (PA) 
approximation applied to the HRM 


The PA approximation of short-time diffusion proper¬ 
ties is based on the cluster expansion of the 7V-particle 
translational mobility tensors of colloidal spheres, 


Mij (X) — 1 Sij + (X) /x 0 1(5^. 

+ three-body terms + ... , 


(Al) 


where the three-body and higher-order cluster contribu¬ 
tions are disregarded. Here, yJ 0 = Dq/^bT) is the single¬ 
particle translational mobility coefficient depending on 
the hydrodynamic particle structure, and 1 is the unit 
tensor. The hydrodynamic mobility tensors on the pair¬ 
wise additive level, 


dV (x) = a4 


N 


Sij I 1+ tJn (R,() | (A2) 

1=1^1 


+ (1 — Sij) U>12 (Rij) 


u>ii and u)i 2 describe the hydrodynamic self-interaction 
of a sphere through flow reflections at the second one, 
and cross-interactions of the two particles, respectively. 
The axial symmetry of the two-sphere problem allows for 
splitting these tensors into longitudinal and transversal 
components, 

<*>ij = Xij{r )rr + j/y (r) [1 - ??] . (A3) 

The transversal and longitudinal mobilities x,j(r) and 
yij(r), with i,j £ { 1 , 2 }, can be calculated recursively 
in the form of a power series in the reduced inverse pair 
distance ah/r, combined with known near-contact lubri¬ 
cation expressions 0 , 0 Ha- 

Insertion of Eq. (IA2I) into Eq. (12U) leads to the PA 
approximation expressions for Dg and Hd(q) [42| . We 
present these here in a form suitable for HRM particles 
where one distinguishes between the hydrodynamic par¬ 
ticle diameter, Uh = 2 ah, and the direct interaction di¬ 
ameter a. For rigid spherical particles, cr = 2a is the 
hard-core diameter, while for mechanically soft particles, 
a is identified with a characteristic soft diameter a s such 
as the one in the Hertz potential. 

Introducing the reduced length x = r/crh and 
wavenumber y = qcr h, the self-part contribution to H(q) 
is 

pOO 

= 1 + 8 y 3 </> / dxx 2 g( r yx) (A4) 

PA J 0 

• [xn (x) + 2y n (x) - 3] . 
The distinct part contribution reads 

(A5) 

"y)' 


J l (xy) . 

■ - J2 ( xy ) , 

L x v J 

where j n is the spherical Bessel function of order n, and 
h = g — 1 is the total correlation function. The over¬ 
lines in xi2(x) and xi 2 (a;) indicate that the respective 
far-field parts up to third order in \/x have been sub¬ 
tracted off. The argument yx in g and h, with 7 = cr^/cr, 
is a reminder that the RDF is commonly calculated as a 
function of r/cr such as in the VW-PY solution for hard 
spheres. 

According to Eq. (1 ATI) , the first-order virial coefficient 
of Dg is given by 

pOO 

A t ( 7 ) = 8 y 3 / dx x 2 exp [—j3V (x)] 

J 0 

• [xn (x) + 2yn (x) - 3] . 


H ( 2 /, 7 ) I PA = 7 0 

poo 

+ 18717 ) / dx xh (yx) 


. ( s 31 (xy) . 2 h(x 
Jo (xy) -+ 7 


xy 


6 x : 


poo 

+ 24y 3 < 7 ) / dx x 2 g (yx) y 12 (x) j 0 (xy) 

Jo 

pOO 

+ 24y 3 <7) / dx x 2 g (yx) [x i2 (x) - y 12 (x)] 
Jo 


Ds (0,7) 
oi ( 7 ) 


with R, j = Rj — Rj , are fully accounted for including the 
near-contact lubrication terms. The two-particle tensors 


(A 6 ) 
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The PA approximation expression for the high- 
frequency limiting viscosity, generalized to the HRM 
reads, 


Voo (</>, 7 ) 


Vo 


PA 


= l + ^[ 1 + 7 3 


(A7) 

/•OO 

60 (y 3 ^) 2 / dxx 2 g ('yx) J (x) , 
Jo 


with the rapidly decaying two-sphere shear mobility func¬ 
tion, J(x), accounting for the two-body His. It decays 
asymptotically as J(x) ~ (15/128)a; -6 . We employ an 
accurate numerical table for J(x) based on rec ursio n ex¬ 
pressions and the lubrication analysis given in Il05l |. 

Finally, we note that in the spherical annulus model 
of hydrodynamically structured hard spheres, the lower 
integration boundary of all integrals in the present 
appendix is equal to I/7. 


Appendix B: Self-part corrected Beenakker-Mazur 
method applied to the HRM 

Except for 7700 , the applicability range of the PA ap¬ 
proximation is restricted to lower (j) values where it be¬ 
comes accurate. In contrast, the ^ 7 -scheme for H(q) and 
Voo by Beenakker and Mazur [70, [Til] is applicable also to 
concentrated systems. We discuss here a standard ver¬ 
sion of this scheme, generalized to the HRM, where S(q) 
is the only required input . On ce the self-part Ds of H{q) 
is suitably corrected [42], Il 06 | , the <5q method provides a 
decent description of short-time transport properties, for 
neutral and charge-stabilized particle systems alike (see 
[Hj]). The <5y method results for 7oo and H(q) reveal 
inaccuracies at all concentrations which can be partially 
attributed to its approximate treatment of the His. This 
is underlined in recent work by Makuch and Cichocki 
[ 72 ] where the approximation steps in the derivation of 
the 67 scheme have been reduced. The fact that their 
revised version of the <5y scheme with improved hydrody¬ 
namic mobility tensors does not significantly improve the 
agreement with simulation data for no-slip hard spheres 
points to a fortuitous cancellation of errors introduced in 
the approximate derivational steps of the original (non 
self-part corrected) Beenakker-Mazur method. 

A simple yet significant improvement over the original 
6"/ scheme preserving its analytical simplicity is obtained 
from using this scheme for the distinct part Hd(q) only, 
where i t giv es overall good results. This was shown in 
Fla . 1671 Il 06 | both for neutral and charge-stabilized parti¬ 


cle systems. Regarding the self-part, Ds, of H(q), accu¬ 
rate expressions can be used instead such as the scaling 
relation in Eq. (BH1) for particles with hard-core interac¬ 
tions, or the PA approximation expression in Eq. (IA4I) 
for lower-concentrated charge-stabilized systems. This 
hybrid procedure is referred to as the self-part corrected 
<5y scheme. 

The ^7 scheme expression for Hd(q) by Beenakker and 
Mazur 70.. 'fH, generalized to the HRM, reads 



with y = qah■ The function S l0 (x) consists of an in¬ 
finite sum of wavenumber-dependent contributions de¬ 
pending on (j)}, as well as on the inter-related scalar co¬ 
efficients with n £ {0,1,2, •■•}. Explicit expres¬ 
sions for S l0 (x) and 7 ^ are given in We have 

calculated the coefficients in an iterative procedure 
up to n = 10 , using a fine grid of volume fractions in 
[0.01 — 0.5] of grid size A (f> = 0.01. This has resulted in 
an improved accuracy as compared to the original work 
by Beenakker and Mazur. However, the differences in 
Hd(q) and 7700 are quite small, i.e. there is no more than 
a 3% difference. 

The £7 scheme expression for 77 ^ by Beenakker . 69} . 
adapted to the HRM, is given by 


Voo ((/>, 7 ) 


Vo 

where 


S-y A(°) {(p, 7) + A( 2 ) ((j>, 7) ’ 


A {0) ( 0 , 7 ) = 


1 + -7 3 </>7q 2) 


A (2) (</>, 7) = ^ 7 3 ^ 7q 2) A (0) (</>, 7) 


1 2 


d y 


ii 2 Gv) 


(B2) 


0 ' I+7 3 J)S 10 (\y, y 3 <(>) L \7 


(B3) 
(B4) 
- 1 


Here, 7 q 2 ^ = 7 q /ti with n denoting the particle number 
density. The argument 77/7 in the static structure factor 
is used as a reminder that it is usually calculated as a 
function of qa. 
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